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I. INTHODDCTION 



A. 07ERVIEW 

Briefly stated, direct time-domain techniques for 
transient radiation and scattering are broadly concerned 
with the response of some object of interest, either 
equipment or a target, to non-monochromatic alectromangetic 
sources. The analysis is done in the time-domain, with time 
treated as an explicit independant variable in. most methods. 

In the past, frequency domain analysis has been the 
primary means of analyzing non-monochromatic SK phenomena. 
Time-domain solutions in closed-form were almost impossible 
to obtain, except in certain very simple cases. The amount 
of computational work required was correspondingly 
significant. No practical aides were available to assist in 
performing the operations prior to the development of fast 
computer systems. 

Another limitation was the requirement that the 
time-domain response be as nearly an impulse as possible and 
thus possess a broad band of frequencies. Hardware was 
required to produce very short pulse waveforms. 
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The situation changed radically with the widescale 
introduction of computers capable of performing rapid 
computations involving large numbers of iterations in 

relatively short periods. This, coupled with the rapid 
maturing of short-pulse hardware technology, soon made 
possible the application of time-domain techniques to a wide 
variety of e letromagnet ic problems. 

.filler and Landt [Kef. 1] list some advantages of 
direct time-domain solutions over frequency domain 

treatments of transient problems. They are: 

1. Greater solution efficiency for many types of EM 
problems. 

2. More convenient handling of non-linearities. 

3. Improved physical insight and interpretablility of EM 
phenomenon. 

4. Availability of wide-band information from a single 
calculation. 

5. Opportunity to isolate interactions, using time range 
gating (e.g., pulse reflections from distant objects 
not of interest, etc.,). 

The major disadvantage is the increased complexity of 
the computer implementation of a time-domain program as 
opposed to one involving frequency domain techniques. There 
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is also a corresponding difficulty in the software 
development and its practical use. Finally, the 
computational time may be significantly greater than that 
required of other techniques. 

The applications for time-domain techniques are diverse. 
Transient analysis can be used to directly evaluate the 
transient characteristics of new antenna designs, 
determining such parameters as driving point currents, 
radiated far fields and near fields, and assisting in 
determining antenna gain. These methods may also be used in 
determining the degree of coupling due to electromagnetic 
pulses (BMP) generated because of lightening during a 
nuclear detonation. Finally, time-domain techniques may be 
directly applied in the solution of the inverse scattering 
problem involving the receiving and scattering properties of 
a large category of target types. 

As time-domain transient techniques are further refined, 
it can be expected that a greater number of uses will be 
identified. 

B. aE7IESr OF TSAHSIEHT ES ANALYSIS TECHNIQOES 

There are three broad categories of techniques 
identifiable in the literature involving transient analysis. 
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They are direct (t ime- dooiain) techniques, transform 

(frequency domain) techniques, and hybrid 

techniques. [Ref. 23 This report is only interested in 
direct techniques and will not review the others. 

1 . Physical Optics Inverse Scattering 

The earliest techniques developed involved physical 
optics inverse scattering. Kennaugh and Cosgriff [Ref. 3] 
in 1958 proposed that the impulse response of a scattering, 
conducting body is simply the second derivative of the 

projected area function of the body if physical optics 
currants are postulated on the surface. If a target is 
illuminated with an incident ramp plane wave, then the 
projected area function of the target is directly 

proportional to the far-scattered ramp response. This is 
only an approximate relationship, however, since it is based 
on the assumption of physical optics currents. 

A more direct and precise ramp response technique 
was suggested by Kennaugh and Moffat [Ref. 4] in 1965. 
They predicted satisfactory approximate ramp response 
signatures could be obtained with narrower bandwidth 

interrogating signals than that required for the impulse 
response. 
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The above techniques vere generally applicable for 
targets which were rotationally symmetric and for whom the 
direction of the electromagnetic field incidence was axial. 
In 1976, Young [Ref. 5] demonstrated that a volume estimate 
can be obtained from the ramp response waveform. He further 
extended the ramp response technique to arbitrarily shaped 
targets with the incident field being other than axial. 

2 . ”Sxact” Inverse Scattering Method 

Another method which improves upon the physical 
optics inverse scattering theory is the '’exact*' inverse 
scattering method first suggested by Bennett et. al., 
[Bef. 6] in 1974, and further refined in 1977. [Sef. 7] 
The inverse scattering problem is solved using an inversion 
of a space-time integral. The method talces into account the 
backscatter response of the target, the direct effect of the 
incident electromagnetic field, and the correction currents 
flowing on the target surface. Theoretically, by completely 
describing the interactions of the rwo components of the 
surface currents, a close approximation of the target shape 
is possible. The actual implementation is limited to 
"perfectly" conducting scatterers. It is this method which 
has been selected for use in the Naval Postgraduate School 
Time Domain Radar Laboratory (TDEL) . 
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3, Other Methods 



Other methods that have been developed, although to 
a less satisfactory degree of accuracy than the inverse 
scattering methods, include radar target identification 
methods. Only a finite number of solutions (targets shapes) 
are allowed as opposed to the nearly infinite number 
possible using inverse scattering. The models of targets 
are stored for parameter comparison with actual returns. 
Certain parametric aspects of a wideband illuminating 
signal, such as natural resonance frequencies and 
polarization of the return, are used. Other similar methods 
use amplitude and phase differences or special properties of 
Rayleigh region scattering. All these methods are limited 
in practical application. 

Still another important and powerful method 
considers multifrequency target illumination using existing 
radars to simultaneously illuminate targets with 7HF, OHF, L 
and S band radar frequencies with circular polarization. 
[Ref, 8] Circular polarization eliminates any polarization 
angle dependence the signal might have. This approach 
allows for the inclusion of all radar system parameters 
(such as antenna gain, noise level, range, etc.) in the 



identification 



algorithm 



The 



technique 



further 



may be 



demonstrates that the target shape information 
obtained with a high degree of accuracy from radar systems 
correctly tailored to the problem of general target 
identification. 

C. REVIEW 3F TDHL PROGRESS iT UPS 

Three requirements were initially identified. First, a 
method had to be chosen from the techniques listed above. 
Secondly, the construction of a physical laboratory to 
perform the required measurements in support of the chosen 
technique had to be designed and built. Thirdly, an 
alogorithm had to be developed employing the selected 
technique for laboratory use. 

As stated above, the "exact" inverse scattering method 
suggested by Bennett was the scheme chosen. To support 
this, work began in January 1980 where evaluation of various 
designs to be used in constructing a Time Domain Radar 
Laboratory at the Naval Postgraduate School was carried out 
by Capt. C. W. Hammond, USMC, as a thesis project. Work was 
essentially completed with the erection of the laboratory 
and interfacing of hardware in September, 1980. [Ref. 9] 

In support of the TDRL, LCDR Heir Morag, Isreali Navy, 
began, as a thesis project, the development of an "exact" 
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1980. 



The 



inverse scattering algorithm in September, 
algorithm was completed in March 1931 [Sef. 10] However, 
the program was developed in FORTRAN, a language not used by 
the TDRL computer system. A translation to BASIC and 
adaptation of the code was necessary to reduce memory 
storage requirements and to ensure computational times were 
of reasonable lengths. 

The final step in the establishment of the TDRL was the 
validation of the laboratory in confirming the work of 
Hammond and Morag. This is one of the principal aims of 
this project. 

D. THESIS OBJECTIVES 

The objectives of this thesis are to; 

1. Develop a working software program which will acquire 
and prepare transient response data from various 
targets for input to the time-domain integral 
equation. 

2. To show that through signal processing, the effects 
of the receiving antenna can be neglected, and the 
response to arbitrary specified exciting waveforms 
can be constructed, or "synthesized", using PFT 
techn igues. 
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3. Validate the TDRL. 



The main effort has been directed toward the development 
of a user oriented computer program which will allow the 
determination of the second objective and prepare the 
transient domain data for input into the numerical 
time-domain integral equation. The shape of the object 
being observed is displayed both before and after correction 
factors are applied. The shape of the object was developed 
from the ramp response of the impulse function applicable to 
the target. The program was also made compatible with the 
input of data to an alternate natural resonance frequency 
target identification technique, using Prony's method, which 
is being developed by LT Demetrius Papaspiridakos, Hellinic 
Savy, as a thesis project. 

This report first describes the physical structure of 
the TDSL. It then derives the relevant time-domain integral 
equations for both the radiation and inverse scattering 
solution for simple axisymmetrical , perfectly conducting, 
closed surfaces in Chapter III. Next, it discusses and 
develops the numerical solutions to the integral equations. 
In Chapter V, the implementation of the complete computer 
code is presented. And finally, measurements that are 
relevant to the validation and confirmation of previously 

developed concepts are reported. 

24 



This thesis is the third to be produced for this major 
project at the Naval Postgraduate School. It marks a point 
of departure in that essentially all the preliminary work 
which will allow for the effective use of the TDRL has been 
performed with validation. The laboratory is, for the most 
part, ready for applied research involving transient 
time-domain analysis. 
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II. PHYSICAL DESCRIPTION OP TDRL 



A. INTHODOCTION 

As stated in Chapter I., the primary thrust of this work 
is to provide a practical working alogorithm to serve as the 
principal link between the acguisition of transient data for 
targets situated on the TDRL scattering range, and the use 
of that data in the inverse scattering algorithm. 

The work performed on this project is a direct extension 
of the preliminary work performed by Hammond in the 
construction of the Haval Postgraduate School Time Domain 
Radar Laboratory (TDRL) , and by Morag in the development of 
the numerical integral equation algorithm for the solution 
of the inverse scattering problem involving axisymmetric 
ta rgets. 

The results obtained by Hammond and Morag are 
essentially complementary. However, the independant 
approaches each took in thair work, particularly in the 
selection of unrelated processing machines for data 
processing, made necessary certain practical modifications 
to the TDRL. It is germane to review the basic work of 
Hammond and discuss those modifications made. Such a review 
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will also allow for the setting of a frame of referance to 
aid in the explanation of the linking algorithm developed 
and results obtained in this project. 

3. PHYSICAL DESCRIPTION OF TDRL 

The physical set-up of the TDRL is shown in Figure 2. 1. 
As can be seen, the TDRL can be broken functionally into 
three groups: 

1. The Impulse Generator Source Group 

2. The Imaging Plane Group 

3. The Sampling Receiver and Signal Processing Group 
Basically, the TDRL is an open system designed to determine 
the transfer function of a target from its backscattered 
radiation. Scattering measurements are bi-static. A 
description of each functional group follows. 

1 . Source Group 

This group consists of a single item, a Video 
Communications OHF Impulse Generator, Model 1000. This 
generator is the source of the transmitted pulse and also 
provides the triggering signal which synchronizes the Signal 
Processing Group to the transmitted signal. 

The Impulse Generator is specified to provide a 
fixed amplitude 45- volt pulse with a maximum rise time of 
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Figure 2.1. Tiae Domain Radar Laboratory 
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Figure 2.3. Source Pulse fro« Generator 

370 picoseconds. Figure 2.3 is a typical waveform as 
measured on a Tektronix S-6 sampler Head by the 
manufacturer . 

The output of the Impulse Generator as measured at 
the TDRL as measured by the Operating System is given in 
Figure 2.4, and Figure 2.5 as photographed from the DPO 
directly. To set full waveform presentation on the measuring 
oscilloscope, 30dB attenuation was applied to the generator 
output. Taking the attenuation into consideration, the 
measured peak value of voltage (1.4604 volts) gives the 
actual output of the generator to be 46.2 volts. This is 
about 2 % better than rated. The rise time of the output is 
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TARGET; RUN DATE DIST ANT TGT REMARKS 

I 6-2S-81 N/A H0 7?7 GENERATOR HAUEFORH 

MAXIMUM PEAK UALUE <UOLTS) 1.1284 

MINIMUM PEAK UALUE <UOLTS) -0.3320 

RMS UALUE (UOLTS) 0.4214 

MEAN UALUE <UOLTS> 0.0808 

HUMBER OF UAUEFORMS AUERAGED > 1 OPTIMIZATION UALUE - 8 




Figure 2.4. Source Wavefora Measured at TDRL 

about 250 picoseconds, also better than specified by the 
manufacturer, by about 48%. 

Some notable features of the generator used in TDRL 
measurements are a variable trigger repetition rate from 10 
Hz to 1 MHz and an adjustable delay from 50 to 150 
nanoseconds. All outputs are matched into a 50-ohm load. 
The generator allows rapid acquisition of highly coherent, 
repetitive waveforms. This is an important feature when 
attempting to maximize signal-to-noise ratio through signal 
averaging. 
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Figare 2.5. DPO aavefocm — OHP Inpalse Generator 

The Video Communications Impulse Generator replaces 
a Tektronix Type 109 impulse generator used by Hammond. The 
completely solid-state electronics of the Mark 1000 has 
proven more reliable in producing a relatively noiseless, 
pulse-to-pulse coherent signal than did the reed-switched 
Type 109. The Type 109 had a fixed pulse repetition rate at 
about 720 pulses per second. This is a very limiting factor 
when attempting to noise average a signal. 

The impulse source is connected to the transmitting 
antenna via a 3. 2 meter RG-8a/0 coaxial cable. This cable 
serves both as a conductor of the source pulse and a short 
delay element to the transraitred signal. 
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2. Imaqipa Plane Group 

The Imaging Plane Group consists of four elements 
the transmitting antenna, the receiving antenna, the 
target, and the imaging plane, 

a. Transmitting Antenna 

The transmitting antenna is a 6.4-meter endfed 
monoDole. It extends vertically from the ground plane and 
produces a vertically polarized EM field. 

b. Receiving Antenna — Transverse Electromagnetic 

Horn (TEM) 

Hammond explored a variety of receiving antenna 
types to be used in the TDRL. The most useful is the TEM 
horn antenna in that it has the following advantages over 
the other antennas considered: 

1. Better directivity. 

2. Less Attenuation, 

3. Mobility on the image plane. 

4. A readily variable characteristic impedance. 

5. ^ Relatively frequency independant over a wideband. 

The TEM Horn used in the TDHL is physically an 
open sided pyrimidal structure when deployed. It is a 
single piece of flat copper, about 1mm thick and triangular 
in shape. Figure 2.5 gives the relevant dimensions. The 
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432 mm 






Figure 2.6. TEa Horn Dimensions 

hole at the apex of the antenna allows for connecting to the 
receiving cable. 

The receiving horn antenna is mounted in a three 
dimensional trapezoidal support mount with quarter-inch 
plexiglass walls on the top, bottom, and sides. The front 
and rear are open to allow easy insertion of the copper 
plate. The dimensions and shape of the mount are given in 
Figure 2.7 All dimensions are outside wall measurements. 

The flat copper antenna is slipped into the 
front of the plexiglass mount and protrudes through the 
back. A coaxial cable is connected to the back with an rf 
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C. FRONT 



Figure 2.7, TEB Horn Support Hount Diaensions 



end-connector and hand-tightened screw. See Figure 2 
a diagram of the assembly of the connection point 
antenna is physically supported within the mount by • 



13 for 
The 
rooves 
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cut into and running along the length of each side. The 
grooves serve to fix the flare angle of the antenna at 
approximately 24. 10. 

The dimensions of the mount were chosen for 
reasons of practicality, being the most appropriate for the 
material available and providing a reasonable flare angle 
and antenna characteristic impedance for signal reception. 
The mount also provides stabilty to the Horn under various 
weather conditions. Thus, by placing the horn in the 
plexiglass support mount, physical and electrical 
characteristics are fixed, assisting in correlating run to 
run measurements. The TBH Horn as it would be mounted on 
the imaging plane is drawn in Figure 2.8. The support mount 
is not shown for clarity. 

The addition of the plexiglass mount should have 
negligible effect on the received signal. Figure 2.9 shows 
the response of the supported horn antenna to a pulse 
generated by the source. This compares favorably with the 

V 

pulses obtained by Hammond without the plexiglass 
mount . [ Bef . 11] 

As noted in Figure 2.8, the electric field at 
the aperture of the horn, assuming a uniform TEM wave being 
excited between the horn and the ground plane, is given by: 
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(a) Front View 




(b) Side View 



Figure 2.8. TEM Horn Heceiving Antenna 
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where: 7^ = amplitude of the incident 7oltage at 

the feed point. 

h = height of aperture 

k - wave number = 2TT/^ 

Tq = distance of observer from the 
reference point. 

It can then be easily shown that the receiving 
antenna impulse response in the far field is; [Bef. 12] 

S^(t) = - c WZ^ 5 (t - r^/c) (2.2) 

where: & = impulse 

= load impedance 

and the impulse response of the transmitting antenna is: 

^ 60 e W 

S^t) = - ^ 5 (t-r^/c) 

where: r = distance from the aperture to an 

observer on the image plane in the 
boresight direction. 

Thus the transmitting impulse response is simply 
the time derivative of the receiving impulse response. The 
impulse response is a delayed impulse, delayed in time by 
r^/c, and multiplied by a negative constant which is 
dependant on the width of the antenna aperture, w, and the 
load impedance Z,_. 
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1 
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e.sa i.ee i.sa a.aa z.sa s.aa 3.sa 4.aa 4.sa s.ea 

Figure 2.9. Typical TEM Horn Received Waveform 

Generally# the following limitations are 
applicable to the TEM Horn during most TDRL 

measurements; [Ref. 13] 

1. The antenna is not in the far field of the 
transmitting antenna. 

2. The flare angle of the TEM Horn is not wide enough to 

receive energy from all positions along the 

transmitting antenna. 

3. Wavelengths at high frequencies may not propagate in 
a true TEM mods# resulting in attenuation and 
distortion. 
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The net result is a deviation 



from the ideal 



results of the impulse responses as stated above. When 
Figure 2.9 is compared with Figure 2.4, distortion of the 
pulse as received by the antenna is noted. However, for 
TDRL purposes, the reproduced waveform is still reasonably 
faithful to the original source pulse. 

c. Receiving Antenna -- Monopoles 

Monopole antennas are also used on the TDRL 
scattering range. One of eleven monopoles, ranging in 
length from 85 millimeters to 285 millimeters, in 20 
millimeter lengths, can be chosen. Each monopole is 
constructed of brass rod and is threaded on one end to allow 
for ease of connection to the receiving elements. 

The monopoles may be considered center-fed 
dipoles with their images in the ground plane. In effect, 
they are oscillating electrical dipoles in free space. 
Oscillation at frequencies related to the 

quarter-wavelengths of the monopoles is to be expected. 

Figures 2.10, 2.11 and 2.12 confirm the above. 

The effects of radiation damping, dispersion and reflection 
due to the impedance discontinuity of the antenna are 
clearly demonstrated. From Figure 2.4 the pulse width of 
the source pulse is seen to be about 1.2 nanoseconds. This 

corresponds to a wavelength of about 360 millimeters. 
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Figure 2,10. 85aa Honopole Pulse Response 





Figure 


2.10 is the short-pulse 


response of 


the 


85mm monopole. 


This antenna 


is almost exactly 


one 


quarter- wave length 


in height in 


relation 


to the source 


pulsewidth. 


It 


correspondingly 


produces 


a nearly 


pure 


sin usoidal 


variati 


on, damped in 


time. 


as predicted 


by 



classical antenna theory. 

Figures 2.11 and 2.12 are the observed responses 
for the 135 lillimeter and 295 millimeter monopoles. As can 
be seen, the antennas tend to resonate at their 
quarter-wavelength periods of 2.5 and 3.8 nanoseconds 



41 



TARGET; RUN DATE DIST fil 
3 6-29-81 1.27M 7 

HAXIMUI1 PEAK UALUE (UOLTS) 

MINIMUM PEAK UALUE (UOLTS) 

RMS UALUE (UOLTS) 

MEAN UALUE (UOLTS) 

NUMBER OF WAUEFORMS AUERAGED » 1 



TCT REMARKS 

??7 ACQUIRE 18SMM ANTENNA RESPONSE 
9.0859 
-0.0907 
0.0329 
9.0000 

OPTIMIZATION UALUE ■ 0 




r aspectively . Although the amplitudes of these two 
aonopoles do not differ greatly (less than .08% between 
their maximum peak values) , the longer antenna tends to 
distort the incident wave to a greater extent. This is most 
readily observed when comparing the root-mean-sguare 
voltages of the two waveforms. The 185 millimeter antenna 
provides more effective power to the receiving circuits than 
does either of the other two monopoles. The 285 millimeter 
antenna dissipates a considerable amount of energy in the 
long period oscillations; the 85 millimeter antenna is too 



small to acquire the full power of the traasmitting antenna. 
There may then be a "best” monopole length for a particular 
source pulse. what this length might be and whether the 
monopoles are selective of targets according to their 
dimensions will be discussed more fully in Chapter 7. 

TARGET: RUN DATE DIST AHT TCT REMARKS 

2 S-29-et 1.27M 12 ??7 ACQUIRE SSSHM ANTENNA RESPONSE 

MAXINUN PEAK UALUE <U0LTS> 9.0661 

MINIMUM PEAK UALUE (UOLTS) -0.0934 

RMS UALUE (U0LTS> 9.0289 

MEAN UALUE (UOLTS) 0.0000 

HUMBER OF HAUEFORMS AUERACED ■ 1 OPTIMIZATION UALUE ■ 0 




Figure 2.12. 285 mb Bonopole Pulse Response 



All receiving antennas are coupled to the Signal 
Processing Group through a coaxial end-coupling fitted 
through the imaging plane. The connection point is 



diagrammed in cut-away in Figure 2.13. 



The antenna is 




INTESIOR TIME DCMAIi'i 
MEASUREIffiNT LABORATORY 



Figure 2.13. 



Beceiving Autenna Connection Point 
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coupled as close to the ground plane as possible. This is 
particularly important for the TEM Horn to prevent 
undesirable oscillations. As drawn^ the threaded portion of 
the monopole (or hand tightened screw of the TEU Horn) is 
joined to a connector pin by butting the surfaces of the 
antenna and pin together. A coaxial cable is then connected 
to the end connector and the rf signal is transferred to the 
receiving circuits. 

This arrangement has worked well. No 

reflections or excessive attenuation of the incoming signal 
due to the joint were observed. Oscillations on the TEM 
Horn were minimal. 

Table I lists the antennas currently available 
to the TDRL user. The significant dimensions and pertinent 
characteristics are reviewed, 
d. Targets 

The targets used on the TDRL scattering range 
are of various dimensions. Currently, eleven are used in 
imaging experiments. They are all axisymmetric, closed 
surfaces with edges and vertices. The targets may be broken 
into three catergories according to general shape: 

1. Half-cylinders 

2. Half- spheres 
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TABLE I 



TDSL Antennas 



ANTENNA 

NOHBER 


DESCRIPTION 


DIHENSIONS 


COMMENTS 


1 


TEH Horn 


Figure 2.6 


Broadband 


2 


Dipole 


85mm 


Oscillating Dipole 


3 


Dipole 


1 05mm 




4 


Dipole 


125mm 




5 


Dipole 


145mm 




6 


Dipole 


165mm 




7 


Dipole 


185mm 




8 


Dipole 


205mm 




9 


Dipole 


225mm 




10 


Dipole 


245mm 




1 1 


Dipole 


265mm 




12 


Dipole 


285mm 





3. Half- cones 

Table II lists the targets and provides general information 
about each. 



There are six half-cylinders used as targets. 
They are all machined from aluminum stock. The largest is 
12-inches long with a 3-inch radius. The smallest is 
3.06-inches long and 1.5-inches in radius. 

A single half-cone is also used as a target. 
The cone is machined from the same stock as the 

half-cylinders. 
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Four spheres, varying in radius from 5-inches to 
1.5-inches are the four remaining targets. The spheres are 
basically styrofoam balls cut in half. Each half is covered 
with heavy-duty aluminum foil to provide the proper 
reflection of transient EM. 

Each of the targets, when placed on the ground 
plane, appears to the incident pulse to be mirrored. 

Therefore, the receiving antenna sees not the half-shapes, 
but the full cylinder, cone cr sphere. 

No particular problems were noted when using the 
targets except for the need to anchor the spheres with 
weights during inclement weather. Otherwise, they have a 
tendency to move on the plane due to the wind, 
e. Image Plane 

The final element of the Image Plane Group is 
the image plane itself. It is a square structure of welded 
Reynolds Type AN-190 aluminum sheet. It measures about 11 
meters on a side for an area of about 121 square meters. It 
is four millimeters thick. 

Based on the above, the range clear time, which 
establishes the low frequency cut-off of the measurement, is 
greater than 36 nanoseconds. The clear time establishes 
the first arrival at the measurement point of a reflection 

from objects surrounding the plane and the end of the plane. 
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TABLE II 



TDRL Targets 



TARGET 


1 

TYPE 


RADIUS 


LENGTH 


KATERIAL 


1 


Cylinder 


3» 


12" 


Aluminum Stock 


2 


Cylinder 


1.5" 


12" 


Aluminum Stock 


3 


Cylinder 


3'’ 


6" 


Aluminum Stock 


4 


Cylinder 


1.5" 


6.0 6” 


Aluminum Stock 


5 


Cylinder 


3" 


3.06" 


Aluminum Stock 


6 


Cylinder 


1.5" 


o 

• 

1 


Aluminum Stock 


7 


Sphere 


5" 




A1 and Styrofoam 


8 


Sphere 


4" 




Al and Styrofoam 


9 


Sphere 


3" 




Al and Styrofoam 


10 


Sphere 


1.5" 




Al and Styrofoam 


11 

1 


Cone 

j 


3.13” 
1 


5.94" 

i 1 


Aluminum Stock 
1 



The main purpose of the image plane is to 
isolate the targets and the antennas from the underlying 
structure of the TDRL. Instrumentation cables are 

introduced from under the ground plane and are thus kept 
from interfering with the measurement. The plane also 
serves as the primary support structure for both the 
transmitting and receiving antenna. 

No significant problems were encountered with 
the image plane. However, some disadvantages of this plane 
arrangement are that only objects having a symmetry plane 
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can be measured, and polarization and incident angles are 
limited. But, for the work intended, the image plane serves 
its purpose well. 

3. Sailing Heceiver and Signal Processing Grou p 

The Signal Processing Group consists of three 
elements: 

1. Digital Processing Oscilloscope (DPO) 

2. Mini-Computer 

3. Hard Copier 

The Signal Processing circuits are connected from 
the output of the receiving antenna to the input of the DPO 
via 2.58 meters of RG-8A/0 coaxial cable. The cable 

provides some attenuation and delay of the source signal 
provided by the impulse generator. Additionally, a trigger 
signal is sent to the DPO from the pulse generator via .54 
meters of RG-53C/0 cable. 

The lengths of the signal and trigger cable are 
chosen for convenience and for minimum attenuation. They 
may be varied for special purpose applications such as to 
accommodate a different trigger source by providing a 
necessary delay. Care must be taken not to introduce too 
much delay or the received signal will not be windowed on 
the CRT display of the DPO. The maximum practical length of 
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cabling, including that to the transmitting antenna from the 
pulse generator, should be less than 20 meters for the 
current equipment configuration of the TDRL. 

a. Digital Processing Oscilloscope 

The Tektronix Digital Processing Oscilloscope is 
a computer compatible, analog to digital system. It is 
comprised of two elemental units; a Tektronix 7704A General 
Purpose Oscilloscope System, and a P7001 Processor. 
Modification of the 7704A Oscilloscope by the addition of 
various plug-in modules allows for tailoring of the whole 
DPO system to the specific needs of the TDRL. A cursory 
description of the various elements of the DPO, as 
applicable to the TDRL, follows. 

The 7704A Oscilloscope system is composed of two 
parts; the D7704 Display Unit and the A7704 acquisition 
unit. The Display unit provides results of signal 
acquisition, both real time and as stored in the P7001. It 
does this by a visual interface with the DPO in the form of 
an Oscilloscope. Additionally, the CRT readout is carried 
over into the DPO to provide appropriate scaling on the 
displays. 

The A7704 acquisition unit, through its modular 
plug-ins, provides the means to acquire the analog signal. 
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The plug-in units used in the TDRL are the S-6 Sampling Head 
and the S-53 Trigger Recognizer, inserted into a 7S12 
TDR/Sampler. 

The 7S12 provides for measurement of recurring 
fast-rise time signals; i.e., the source pulse and the 
target bacJc-scatter. It features the means of determining 
vertical deflection factors in millivolts; horizontal sweep 
factors in seconds; a HIGH RESOLUTION switch which reduces 
the waveform noise and jitter by signal averaging (10 
averages per sample point) ; a TIHS-DIST ANCE scale which 
allows one-way distance measurements in air dielectrics; and 
a locate switch which increases the tine/division and 
intensifies a portion of the display to locate the time 
window relative to the total waveform. Finally, the 7S12 
provides the means for the inserting and interfacing of the 
plug-in heads into the DPO. 

The S-6 Sampling Head is a 50-ohm loop through 
input sampling unit. It provides the means of inputting the 
signal of interest to the DPO. This unit has a specified 
rise time of 30 picoseconds or less and a bandwidth 
equivalent to dc to 11.5 GHz at 3 dB down. It has a signal 
range of ♦ Ivolt to -1volt (dc to ac peak) or Ivolt 
peak-to-peak . The S-6 establishes the baseline measuring 
capabilities of the TORL. 
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The S-53 Trigger Recognizer permits use of the 
7312 as a general purpose sampler. It produces a stable 
trigger signal from input signals from DC to 1 GHz. It is 
capable of recognizing signals having amplitudes ranging 
from 10 millivolts to 2 volts peak-to-peak into 50 ohms. A 
trigger pulse from the pulse generator is fed into this unit 
at the same rate as the source signal. Trigger to signal 
delay is specified as 15 nanoseconds or less. 

The second element of the DPO, the P7001 
Processor, provides the link from the DPO to the 
mini-computer. This unit receives the analog signal 
acquired from the Acquisition unit and digitizes and stores 
or discards the information, along with appropriate scale 
factors, as selected by the operator. The method used by 
the digitizer is pse udo- random sampling with one sample 
taken every 6.5 microseconds. This is the maximum rate of 
sampling. Any transient longer than 5 milliseconds or 
COHERSNTLT REPETITIVE signal up to the frequency limit 
specified by the Sampling Head that can be displayed on the 
D7704 Display Onit can be stored in internal memory along 
with its scale factors . The information stored can then be 
redisplayed on the CRT and/or sent to the minicomputer for 
processing. There are a maximum of 512 samples taken 
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horzontally (9-bit word) per waveform. Although either X- 
or 7- (TIH2) may be specified for this axis, the TDRL 
primarily employs the time axis. Vertical resolution of the 
Processor is 1024 levels (10-bit word). This axis is used 
to measure the voltage of the back-scattered signal. 

b. Tektronix 4052 Mini-Computer 

The minicomputer used in the TDRL is a Tektronix 
4052. It is a high performance, integrated graphics system 
employing LSI, bi-polar, 16-bit technology. The unit has a 
64k-byte memory with a 300k-byte magnetic cartridge built 
in. Extended BASIC is the high level programming language 
employed to perform processing functions. Although this 
unit has many features worthy of note, the above are the 
salient ones for the TDRL. 

c. Hard Copier 

The Hard Copier used is the Tektronix 4631 Hard 
Copy Unit, It provides permanent, dry copies of the graphic 
and alphanumeric information displayed on the CRT storage 
screen of the 4052 mini-computer. 

All units described for the Signal Processing 
Group were especially selected to perform as an integrated 
system requiring minimum interfacing. The GPIB bus is the 
means of transferring information from one unit to the other 
in the TDRL. 
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III. THEORY 



This chapter is devoted to the general derivation of the 
equations which are applicable to the tiae-domain analysis 
of the impulse and ramp responses of axisymmetric, metallic 
bodies which are illuminated by an axially directed incident 
electromagnetic field. Where sources providing detailed 
derivations are numerous and/or adequate, the derivations in 
this chapter will be general, providing only the major 
results of interest with specific constraints noted. 

The time-domain integral equations which provide 
solutions for the transient electromagnetic problem will be 
considered first. They will be followed by a discussion of 
the measurement equations and procedures employed in the 
TDRL. Finally, various methods used to reduce the noise 
error in acquired signals will be discussed. 

A. TIHE-DOHAIN INTEGRAL EQUATIONS 

One primary problem in transient analysis of target 
impulse responses is to determine the shape of an unknown 
object when the incident field and the response of the 
scatterer to that field are a priori information. Numerous 
time-domain techniques have been developed which perform 
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this function. As noted in Chapter I, the earliest method 
used involved the physical optics response of a body to an 
incident electromagnetic field. The surface currents 
produced an approximate impulse response that was 
essentially the second derivative of the projected area 
function of the scatterer. [Bef. 5] However, this approach 
was limited in its accuracy in that the relationship between 
the impulse response and the derivatives is exact only at 
the leading edge of the scattered field response, a single 
point in time. Interactions between currents on the target 
body which continue to radiate for a significant time after 
the incident field leading edge has propagated further in 
space alters the backscattered response of the target. 
Therefore, for an accurate solution to the problem, these 
additional currents must be of consideration. 

Bennett [Ref. 14] first proposed a direct method which 
'•corrected" the response by computing the currents flowing 
on the scatterer surface. From these currents, the 
scattered field could be calculated "exactly". 

In general, the remainder of this section follows the 
format of the method used by Bennett in the derivation of 
inregral equations for exact solutions. 
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1 t Integral Solution 

The expression for the aagnetic field H at an 
arbitrary point in space (not on the scatterer surface) is 
given by: 









X a 



R 



ds ' 



(3.1) 



where: a(r,t) = total magnetic field at (r,t) 

{?/ t) = incident magnetic field at (r,t) 
J(r*,t) = surface current at 

r = position vector to the 

observation point 

r' = position vector to the 

integration point 
R = |r-'r' 1 




t = time in light-meters (one 
light-meter is the time it takes 
light to travel one meter) 

T = t-R = retarded time 

3y specializing the arbitrary space point to "r, a 
point on the surface of the scatterer, and then applying the 
boundary conditions (i.e., by causuality, incident field is 



56 



zero prior to arrival of the incident pulse at a sample 
point), an integral aquation for the curreut density on the 
surface of the scatterer can be found; 



= 2a^ X H"(r,t) - ^ ^ x {[4 > | A] 



X a^^} ds'" 



(5.2) 

where: "r = position vector to a point located on 

the .surface of the scatterer 
“a^ = unit vector normal to the surface 
The first term on the right-hand side of equation 
(3.2) is the source term. It represents the direct 
influence of the incident field on the current at the 
observation point (^,t) . It is the physical optics 
approximation for the surface current when applied to the 
illuminated side of the scatterer. The integral on the 
right side of equation (3.2) represents the influence of 
currents at other surface points on the current at (?, t). 
Mote that the influence of other currents on the current at 
(^,t) is delayed by R. This allows the surface current 

density equation to be solved by an iterative numerical 
procedure in the time-domain, rather than the familiar 
frequency-domain matrix inversion process. 

Once the surface currents have been found, the 
far-scattered field can be calculated by using; 
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(3.3) 



r 

o 



H^(r,x) 



4it 




-i- '' 

(r ' , t) X a^} ds 



where: r^ = distance to the far-field observer 

- unit vector from the integration point 
to the far-field observer 

With the substitution of the surface current 
expression into the above equation, and the assumption that 
the incident field is a ramp, the result for the backscatter 
direction is: 



■ ST * W 17 ^ 



X a^} ds 



( 3 . 4 ) 



where: r 

O 



^(tjl 



backscatter ramp response of the 
target 

physical optics silhouette area of 
the target 

distance of the far-field observer 
from the origin 




^ |H^| 

, correction currents resulting 
from the incident ramp waveform. 
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This result provides the exact relationship between 
the target response and target geometry. A particular 
unknown sample value at a given point in space and time is 
determined by the exciting field at that same space-time 
point and by the scattered fields from earlier, more distant 
locations. Note that the interactions between unknown 
samples are displaced in time by an amount equal to the time 
required for a field to propagate between the samples at the 
speed of light. The unknown samples can be solved at any 
time step, provided all sample values at earlier times are 
already known. Thus, by determining the unknown samples, 
i.e., the correction currents '3'^, and adding the 
contribution of the target area function, s(t^) , the target 
ramp response can be utilized to recover information 
involving the unknown dimensions. 

2* In verse Scatte rin g Solution 

The targets employed in the TDRL are all 
rotationally symmetric scatterers, similar to the one 
diagramed in Figure 3.1. The general constraints applicable 
to the problem are that the scatterer is symmetric about the 
z-axis, the incident field is axially incident, and the 
far-field is computed in the backscatter direction. 
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The contour of a rotationally symmetric target is 
completely described by the radius vector p (z) which 
varies as a function of z. The projected area function in 
the plane orthogonal to the direction of the propagation of 
the field (z-axis) can be expressed simply as; 

S = TIP^(Z) 

Substituting the projected area function into the 
equation for the backscatter ramp response of the target and 
solving for (z) , an inversion equation for the 
rotationallly symmetric case is obtained: 



p(2) = [2r^H|(?,t) - ^ 4 






a^} (3.5) 



Thus, p (z) is given in terms of the measured ramp response 
(a field measurement) , and in terns of correction currents 
at earlier times, which have been previously computed or are 
known to be zero. 

In the TDRL, an iterative approach is employed using 
estimates of the entire target geometry to solve for the 
contour function. The procedure involves five steps: 

1. Neglecting the surface integral providing the 
correction currents, 
estimate of Q (z) . 



obtain 



the physical optics 



2. Determine the correction currents using the numerical 
solution of equation (3.2) for surface current 
density. 

3. Apply the correction currents to equation (3.5) to 
obtain a next estimate of P (z) . 



4. 



PlCz) = [2r^ h| 



Compare the new value ^ (z) with ^ (z) to see if 



the change is less than some small number, an error 
factor. The normalized mean square error is defined 
as: 






where : 




(t) 



the value to be minimized 
in order to obtain the 
most accurate shape 
description 

the Jcth estimation of the 



backscatt ered field. 

Thus, the difference between the estimated and actual 
backscatt ered fields are directly linked to the 
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difference bet’^een successive contour functions, and, 
therefore, successive shapes of the body. 

5. Repeat the iteration process until a minimum error is 
obtained. 

The above procedure was used successfully by Morag 
[Ref. 15] in the development of an algorithm involving 
theoretical analysis of targets. This algorithm has been 
integrated with the overall TDRL algorithm for inverse 
scattering. 

B. BACKSCATTERED iA?EE0R8 HEASaRESEHTS 

The TDRL algorithm is designed to minimize the effects 
of the impulse responses of the transmitting and receiving 
antennas. Figure 3.2 diagrams the parameters involved. 

As can be seen, there are two measurement paths of 
interest. The first is the Direct Path between the 
transmitting and receiving antenna. The only factors 
affecting the transmitted waveform are the intervening 
distance and the impulse responses of the receiving and 
transmitting antennas to the source pulse. The second path 
is the Augmented Path, from the transmitting antenna to the 
target and back to receiving antenna. Factors effecting the 
transmitted pulse in this path include the antenna impulse 



53 



•iPWiiPri 











IMPULSE 




WIRE 


GENERATOR 


V 


ANTENNA 


(S(w) ) 




(H^(w) ) 



FREE 

SPACE 

DISTANCE D 



(H^( w)) 



FREE 








FREE 


SPACE 




TARGET 




SPACE 


DISTANCE D- 


^ 




^ 


DISTANCE 2D 






'(H^(w)) 






(H^(w) ) 








(H^Cw)) 



DIRECT 

PATH 




RECEIVING 

ANTENNA 

(H^(w)) 



NOTE: Effects of wire transmitting antenna 

reflections at receiving antenna considered 
small — neglected. 



Figure 3.2. 



Signal Paraaeters for TDRL Measurements 





responses, the impulse response of the target and the 
corresponding intervening distances. 

Three waveforms are measured: the Direct Waveform, the 

Incident Waveform and the Augmented Waveform. During the 

I 2 

measurements, h (t) and h (t) , responses due to the delays 
of the intervening distances, are windowed out by using the 
time-range gating feature of the DPO. 

The Direct Waveform is measured first. The transmitting 
and receiving antennas are boresighted. No target is on the 
imaging plane. The measured wave is then the convolved 
responses of the source, and transmitting and receiving 
antenna impulse responses. By using the properties of 

convolution, the frequency domain representation converts 
the convolution process into a product. The frequency 

domain representation cf the Direct Waveform is: 

where: A(e^'^) = Direct Waveform 

S (e"^ ) = source pulse 

= transmitting antenna response 
H*(e^‘^) = receiving antenna response 

The effects of the intervening distance have been time 
gared out. Figure 3.3 is a representative Direct Waveform 
measured in the TDRL. 
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TARGET: RUN DATE DIST ANT TGT REMARKS 

1 8-11-81 1.27M 1 1 CYLINDER ACQUISITION AUERACE 

MAXIMUM PEAK UALUE +289.21 nU 

MINIMUM PEAK UALUE -95.55 hU 

RMS UALUE +90.58 nU 

MEAN UALUE 0.00 nU 

NUMBER OF WAUEFORMS AUERAGED =« 21 OPTIMIZATION UALUE >=0.5 




0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 nS 



Figure 3.3. Typical Direct Savefora--Cylinder 
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The Incident Waveform is measured next. The imaging 
plane is as set-up for the Direct Waveform, but the 
time-gating is further delayed by ddt to eliminate the 
effects of the distance intervening between the receiving 
antenna and the target site. The resulting measured 

waveform is the Direct Waveform delayed by i^t: 

where: 3(6^**^) = Incident Waveform 
Figure 3.4 provides a representative time-domain Incident 
Waveform measured in the TDSL. 

The final waveform measured is the Augmented Waveform. 
It is a convolution of the Direct Waveform with the impulse 
response of the target and summed with the Incident 

Waveform. For this measurement, the transmitting and 

receiving antenna remain boresighted. A target is placed on 
the image plane. The measured frequency domain 

representation is: 

C(e>) = [BCej*^) > A(ej“) 

where: CCe^***) = Augmented Waveform 

H (e^***) = Target Response 

Figure 3.5 provides a representative time-domain Augmented 
Waveform measured in the TDHL for a axisymmetric 
half-cylinder. 
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TARGET: RUN DATE DIST 
1 8-11-81 1.27M 

MAXIMUM PEAK UALUE 

MINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

HUMBER OF WAUEFORMS AUERAGED » 



ANT TGT REMARKS 

1 I CYLINDER ACQUISITION AUERACE 
.. ■»4.40 nU 

.. -3.83 nU 

.. ♦1.92 mU 

. . 0.00 nU 

21 OPTIMIZATION UALUE » 0.5 




0.50 1.00 l.SO 2.00 2.50 3.00 3.58 4.00 4.50 5.00 nS 



Figure 3.4. Typical Incident SJavefor* — Cylinder 
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target: run date dist ant tgt remarks 

1 8-11-91 1.27H 1 1 CYLINDER ACQUISITION AWERAGE 

MAXIMUM PEAK UALUE ♦4.89 nU 

MINIMUM PEAK UALUE -9.76 nU 

RMS UALUE +2.76 nU 

MEAN UALUE 0.08 hU 

NUMBER OF HAUEFORMS AUERAGED » 21 OPTIMIZATION UALUE • 0.5 




a.sa i.aa i.so 2.aa 2.sa 3.ea 3.sa 4.aa 4.sd s.aa ns 



Figure 3.5. Typical ingaeoted Waveform — Cylinder 
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As can be seen from equation (3.6), the target response 
is not independent of the effects of the impulse responses 
of the antennas. To eliminate this dependence, the 
backscattered waveform is determined and then the remaining 
antenna terms are theoretically eliminated by a division 
process. 

The Backscattered Waveform is the difference of the 
Augmented and the Incident Waveforms: 

D(oj“) = C(ej“) - 

= A(e>) (e>) 

Figure 3.6 is the time-domain Backscattered waveform 
resulting from the subtraction of the waveforms in Figures 
3.5 and 3,4 respectively. 

By next dividing the Backscattered Wave by the Direct 
Wave, all antenna influences are eliminated: 

A(eJ“) 

c 

The target impulse response, h (t) , may now be determined by 
deconvolution of the target impulse frequency domain values. 
Figure 3.7 is a representative time-domain target impulse 
response for a cylinder. Noise errors have been minimized 
by an adaptive filtering method to be described shortly. 
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TARGET: RUH PATE OIST ANT TGT REMARKS 

1 8-11-31 1.27M 1 1 CYLINDER ACQUISITION AUERACE 

MAXIMUM PEAK UALUE •*•4.65 nU 

MINIMUM PEAK UALUE -4.44 nU 

RHS UALUE •*•1.71 nU 

MEAN UALUE 0.00 mU 

HUMBER OF WAUEFORMS AUERAGEO « 21 OPTIMIZATION UALUE « 8.5 




0.58 1.88 1.58 2.88 2.58 3.88 3.58 4.88 4.58 5.88 nS 



Figure 3. 6. 



Typical 3ackscattered Wavefora — Cylinder 
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TARCET; RUN DATE DIST ANT TGT REMARKS 

1 8-11-81 1.27M 1 1 TYPICAL CYLINDER RESPONSES 



MAXinUM PEAK UALUE 0.8004 mU 

MINIMUM PEAK UALUE -0.0004 nU 

RMS UALUE 0.0002 nU 



MEAN UALUE 0.OO00 nU 

NUMBER OF HAUEFORMS AUERAGED = 1 OPTIMIZATION UALUE - 10 



nU 



0.0004 



0.0002 



0.0000 



0.0002 




0.S0 1.00 1.S0 2.00 2. 30 3.00 3.50 4.00 4.30 5.00 nS 



Figure 3.7. rypical Target Impulse Response With Noise 
Filtering — Cylinder 
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C. MOISE CONSIDERATONS 

As noted in the previous section, deconvolution of the 
target backscattered impulse response, less the impulse 
responses of the antennas, is a central operation to the 
effective use of the TDRL. The reliability of the acquired 
signal to be a true response of a target to an incident 
field bears directly on the accuracy of the final results as 
determined by solution of the inverse integral equations for 
time-domain analysis. Any errors in the data representing 
the Incident Waveform and Augmented Waveform degrade the 
cancellation of the division process and introduce errors in 
the backscatt ered target response. Also present at the 
input to the DPO is noise due to thermal electron excitation 
in the transmission lines and antenna surfaces as well a 
that due to external atmospheric, cosmic and man-made noise 
sources being received by the horn antenna. When the 
backscattered impulse response is transformed to the 
time- domain, high frequency noise in most cases forms the 
dominant feature of the waveform. An example of a typical 
impulse response of a cylinder is given in Figure 3.8. 
Noise errors have not been removed. Figure 3.9 is the ramp 
response for the same target. 
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target: run date dist ant tgt remarks 

I 9-11-31 1.27H 1 1 TYPICAL CYLINDER RESPONSES 



MAXIMUM PEAK UALUE 0.0324 mU 

MINIMUM PEAK UALUE -0.0314 nU 

RMS UALUE 0.012S nU 



MEAN UALUE 0.O0O0 nU 

NUMBER OF WAUEFORMS AUERAGED » 1 OPTIMIZATION UALUE » 0 



nU 



0.0200 



0.0000 



- 0.0200 



0.50 1.00 1.S0 2.00 2.30 3.00 3.50 4.00 4. SO S.00 nS 



IMPULSE RESPONSE 




Figure 3.8. Impulse Response — Noisy Signal 
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TARGET: RUN DATE DIST ANT 
1 3-11-81 1.27M 1 

MAXIMUM PEAK UALUE 

MINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

NUMBER OP WAUEFORHS AUERAGED » 1 



TGT REMARKS 

1 TYPICAL CYLINDER RESPONSES 
+53. nU-M 
-700.00 hU-M 
+373.03 hU-H 
-285.42 nU-M 

OPTIMIZATION UALUE » 8 




0.58 1.88 1.58 2.88 2.58 3.80 3.58 4.08 4.58 5.08 nS 



Pigare 3.9. Ramp Response — Noisy Signal 
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In the TDRL, typical sources of errors in the Incident 
and Augmented Waveforms are sampling errors and noise in 
signal acquisition, leakage and aliasing errors in the 
Fourier Transf omations , and rounding errors in the computer. 
The error of primary concern is that involving noise in 
signal acquisition and in particular the enhanced 
high-frequency noise resulting from deconvolution of the 
antenna responses to obtain the impulse and ramp responses 
of the target. This noise has been directly traced to the 
structure of the data within the input array, as will be 
noted in Chapter IV. 

Two methods are employed to reduce noise to acceptable 
levels; acquisition averaging and an optimal compensation 
technique. [Ref. 16] 

1 . Ac quisition Av era g in g 

Acquisition averaging is aimed at reducing signal 
acquisition errors and noise by averaging several 
acquisitions taken within a short period of time. During 
the process, a more serious error may be generated, that due 
to signal drift with time. Signal averaging has been most 
useful in ensuring a good target ramp response. It has had 
noticeable effect on making the information content of the 
impulse response more available. Figure 3.10 demonstrates 
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an impulse waveform with 210 acquisitions for the same 
target as for Figure 3.8. Iheoreticaily , the RMS noise 
level should be reduced by the square root of the number of 
acquisitions, N. M = 10 in Figure 3.8, and N = 210 in 
Figure 3.10. The noise reduction should be about 4.5 times. 
Mote that the noise level in Figure 3.10 is indeed about one 
fourth that in' 3.8. 

Figure 3.11 is the corresponding ramp response. 

Signal acquisition only is employed. 

2. Qp t_i ma 1 Compensation Tec hni que 

The optimal compensation technique involves the 
design of a compensator (deconvolution) function operating 
on the convolution output of the Incident and Backscattered 
Waveforms, to produce a relatively noise-free estimated 
Impulse Response for the target. The compensator is applied 
to the convolution product of the target impulse response 
and the Direct Waveform to produce the required 

deconvolution result. The design as used in the TDRL 

involves an iteration process on a single variable with a 
" man-in-the- loop” to determine optimal values of the 
variable. 

Figure 3.12 diagrams the design of a frequency 

domain optimal compensator. The compensation principal in 
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target: run date dist ant tct remarks 

I 8-11-81 1.27M 1 1 CYLINDER ACQUISITION AVERAGE 



MAXIMUM PEAK UALUE 0.0894 nU 

MINIMUM PEAK UALUE -0.0091 mU 

RMS UALUE 0.0033 nU 



MEAN UALUE 0.0000 dU 

NUMBER OF WAVEFORMS AVERAGED * 21 OPTIMIZATION UALUE » 0 



nU 



0.0050 

0.0000 

-0.00SO 



IMPULSE RESPONSE 




0.SO 1.00 1.50 2.00 2.50 3.80 3.50 4.00 4.50 5.09 nS 



Figure 3.10. lapulse Response: No Optinization — Acquistion 
Averaged 
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target: run date dist ant tgt remarks 

4 3-11-81 1.27M 1 1 CYLINDER-ACQUISITION AUC 

MAXIMUM PEAK UALUE ♦82.97 nU-M 

MINIMUM PEAK UALUE -441.33 nU-M 

RMS UALUE +253.42 mU-M 

MEAN UALUE -178.46 nU-H 

HUMBER OF HAUEFORMS AUERAGED * 21 OPTIMIZATION VALUE * 8 




8.S8 i.aa i.sa z.aa z.sa 3.aa 3.sa 4.aa 4.s8 s.aa ns 



Figure 3.11. Raap Response: No Optiaization — Acquistion Av- 
eraged 
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Figure 3.12. Fregueucy Domain Optimal Compensator Design 

deconvolution is to design a transfer function Q (e ) to be 
applied to the backscatter response D(e‘^‘*') to yield an 
estimate of the target impulse response H. (e^***) . That is; 






= D(e>)Q(e>) 



(3.7a) 



* A(ej"')H^(eJ'“)Q(ej“) 



(3-. 7b) 



The main design criteria is to minimize the energy 
in the noise terms that introduce error in the estimate. The 
error energy, B is given by; 

=e • -g” I hJ (e^“) - (.i“) (j.S 

where; frequency band of interest. 

When the error energy is minimized, H^Ce^*^) will 

equal and the impulse responses will be identical, 

except for the addition of highly enhanced random noise. In 

this case, the compensator transfer function is; 
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1 



Hg (e^“) 



(3.9) 



Q(e 




ACeJ^") 



But, this will result in an infinite, unbounded 7alue for 
0(6^*^) as A goes to zero, yielding a noisy form for 

H^(e^*^). This noise can be limited by constraining Q(e ) 
to be a bounded function. Noting that the impulse response 
of the target is bounded by physical constraints, it is 
possible to keep the energy due to the convolution of the 
impulse response and the designed transfer function bounded, 
i. e . : 

Ec = I H(ej“) Q(e^^) dcu (2-103 

The problem now is to minimize E while keeping E 
finite. The total energy in the backscattered wave can be 
defined as: 

E = + XE^ , X _> 0 (3.11a) 

= / [1 (e>) - (ei^) * X i H (e^ ‘") Q (e^ ‘") | ^ ] do) (3.11b) 

O 6 

where; optimization parameter. 
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Thus, by minimizing the total error energy through 

the selection of the optimal value of Xt -he error energy 

in the estimate is also minimized. The result yields an 

optimum compensator for the unknown target impulse 

0 

response h (t) . The result is that the whole system output 
is low-pass filtered with minimal effect on waveform shape. 

Adhering to the design constraints listed above, it 
can be shown that the form of the compensating transfer 
function Q (e^“* ) for a given Direct Waveform A(e^‘^) is: 
[Ref. 17] 



Q(e 







(3.32) 



where: A* (e ) = complex conjugate 

From equation 3.7, the compensation deconvolution 
process yields: 



S _ (3.13) 



Figure 3.13 shows the effect that optimization has 
on the impulse response with no sgnal acquisition averaging. 
Optimization is 0.5. The target impulse response, while 
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- I 



target: ruh date dist ant tgt remarks 

1 3-11-31 1.27M I I typical CYLINDER RESPONSES 



MAXIMUM PEAK UALUE 0.0011 nU 

MINIMUM PEAK UALUE -0.0010 mU 

RMS UALUE 0.0004 nU 



MEAN UALUE 0.0000 nU 

HUMBER OF WAUEFORMS AUERAGED = 1 OPTIMIZATION UALUE *0.5 



nU 



0.0010 

0.0005 

0^0000 

-0.0005 



IMPULSE RESPONSE 




Figure 3.13, 



Impulse Response: .5 Optimization 
No Acquisition Averaging 
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still noisy, is clearly discernable. Figure 3.13 should be 
compared with Figure 3.11. Figure 3.14 shows the effect of 
optimization increased to 10. The noise is almost totally 
eliminated. The trade-off is that the optimization has 
reduced the amplitude of the response and spread it in 
frequency. This is a general observation. However, it 

appears that the degree of effect on the waveform is a 
function of the ideal optimization factor. Chapter V will 
explore this concept more fully. However, it should be 
noted that if the optimization factor is large enough to 
reduce the large scale features of the target then it is too 
large and should be reduced — this is the ’’man in the loop's” 
job. 

The combination of waveform averaging and 
optimization on the impulse response is demonstrated in 
Figures 3,15 and 3.15. With a relatively small amount of 
optimization, the target impulse response is readily 
apparent. See Figure 3.16. Signal acquisition averaging 
also yields an important advantage. Since noise is reduced 
directly by the square root of the number of acquisitions, 
and because no filtering of any sort is employed for 
averaging alone, the detailed features of the actual signal 
are not eliminated as with filtering. The smaller amount of 
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optiaization needed results in minimum effect by the 
compensator on the shape of the output waveform. The 
trade-off is the increased liJclihood of error introduced due 
to signal time drift during the acguistion averaging period. 

Figures 3.17 to 3.21 show the effects of signal 
averaging and optimization on the Ramp Responses. Each 
figure is presented in the same order as were the Impulse 
Responses. Signal averaging is very beneficial in obtaining 
an accurate ramp response. Optimization has a greater 
tendency to effect the wave shape of the ramp response than 
the impulse response of the target. But less optimization 
is required to reduce the error noise. It is generally 
possible to use the ramp response without optimization. The 
main effect of optimization is to spread the frequency 
content of the ramp response. Amplitude effects are minor. 
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target; run date dist ant tgt remarks 

1 8-11-91 1.27N 1 1 TYPICAL CYLINDER RESPONSES 



MAXIMUM PEAK UALUE 8.0QQ4 nU 

MINIMUM PEAK UALUE -0.0Q04 nU 

RMS UALUE 0.0802 tiU 



MEAN UALUE 0.0000 mU 

NUMBER OF HAUEFORMS AUERAGED = 1 OPTIMIZATION UALUE » 10 



nU 

0.0004 
0.0002 
8.0000 
-0.0082 

0.S0 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 nS 



IMPULSE RESPONSE 




Figure 3.14. Impulse Response: 10 Optimization — No Acquisi- 
tion Averagxng 
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TARGET; RUN 0(5TE DIST ANT TGT REMARKS 

I 8-11-81 1.27(1 1 1 cylinder ACQUISITION AVERAGE 



MAXIMUM PEAK VALUE 8.0894 nU 

MINIMUM PEAK VALUE -9.0991 nV 

RMS VALUE 9.9933 nV 



MEAN VALUE 9.0999 nV 

NUMBER OF WAVEFORMS AVERAGED = 21 OPTIMIZATION VALUE » 9 



mV 



9.0959 



8.9009 



0.8959 




9.50 1.00 1.50 2.00 2.50 3.88 3.50 4.00 4.59 5.00 nS 



Figure 3.15. 



Inpulse Response: No Optiaization — 210 Acqui 
sitions Averaged 
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target: run date oist ant tgt remarks 

1 8-11-81 1.27M 1 1 CYLINDER ACQUISITION AUERAGE 



MAXIMUM PEAK UALUE 0.8018 nU 

MINIMUM PEAK UALUE -8.0009 hU 



RMS MALUE 8.0003 nU 

MEAN UALUE 8.0000 nU 

HUMBER OF WAUEFORHS AUERAGED * 21 OPTIMIZATION UALUE « 0.5 



nU 

9.0010 
0.0005 
9.0000 
-0.0005 

0.50 1.98 1.50 2.00 2.50 3.00 3.50 4.00 4.58 5.00 nS 



IMPULSE RESPONSE 




Figure 3.16. Impulse Response; .5 Optimization — 210 Acguis- 
tions Averaged 
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TARGET: RUN DATE OIST ANT 
1 8-11-91 1.2TN 1 

MAXIMUM PEAK UALUE 

MINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

NUMBER OF MAUEFORMS AUERAGED ^ 1 



TCT REMARKS 

I TYPICAL CYLINDER RESPONSES 
+53. GS nU-M 
-TOO. 00 hU-M 
+373.03 nU-M 
-285.42 tiU-H 

OPTIMIZATION UALUE * 0 




0.59 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 nS 



Figure 3.17. Raap Response: Ho Optimization — No Acquisition 
Averaging 
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target: run date dist ant 

1 8-11-81 1.27M 1 

MAXIMUM PEAK UALUE 

MINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

NUMBER OF HAUEFORMS AUERAGED = 1 



TGT REMARKS 

1 TYPICAL CYLINDER RESPONSES 
♦47.77 nU-M 
-847.48 hU-M 
♦367.11 nU-M 
-279.80 mU-M 

OPTIMIZATION UALUE « 0.5 




0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 nS 



Figure 3.18 



Raap Response 

averaging 



5 Optiaization — No acquisition 



target: run date dist ant 

I 8-11-81 1.27M 1 

MAXIHUH PEAK UALUE 

HINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

HUMBER OF UAUEFORHS AUERAGED » I 



tgt remarks 

1 CYLINDER ACQUISITION AUERAGE 
♦69.81 nU-M 
-401.43 nU-M 
♦231.11 mU-M 
-170.58 tiU-M 

OPTIMIZATION UALUE « 10 



nU-H 

50 

0 

-50 

-100 

-150 

-200 

-250 

-300 

-350 

-400 



0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 nS 




Figure 3.19. Raap Response; 10 Optimization — No Acquisition 
Averaging 
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ThRCET: rum date oist ant 

I 9-11-81 1.27N 1 

HAXINUM PEAK UALUE 

HINIMUn PEAK UALUE 

RMS UALUE 

MEAN UALUE 

NUMBER OF WAUEFORMS AUERAGED = 21 



TGT REMARKS 

1 CYLINDER ACQUISITION AUERACE 
+78.58 nU-fl 
-457.95 PiU-M 
+239.59 mU-H 
-223.75 nU-H 

OPTIMIZATION UALUE = 8 




8.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 nS 



Figure 3.20 



Samp Response: 
tions Weragea 



Ho Optimization — 210 Acguisi- 



target: run pate DIST ftNT 

1 3-11-91 1.27M 1 

MAXIMUM PEAK UALUE 

MINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

NUMBER OF HAUEFORMS AUERAGED » 21 



TGT REMARKS 

1 CYLINDER ACQUISITION AUERACE 
+72.25 nU-M 
-465.94 nU-M 
+288.01 fiU-N 
-222.20 rtU-ti 

OPTIMIZATION UALUE • 0.5 



hU-M 

0 

-200 
-400 

0.50 1.00 1.50 2.00 2.50 3.90 3.50 4.00 4.50 5.00 nS 




Figure 3.21. Ra«p Response: .5 Optimization — 210 Acquisi- 
tions Averaged 
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IV. OPERATING SYSTEM COMPUTER IMPLEMENTATION 



A. INTRODUCTION 

A set of software programs called the ''Operating System" 
controls the input, processing and output of information 
obtained in the TDRL. It is the quality and throughness of 
the Operating System that directly determines the quality of 
the results obtained. In this chapter, there will be a 
discussion in detail of the implementation of the Operating 
System. It will begin by discussing the physical structure 
of the the host processor. A brief overview of the 
Operating System will follow, defining the basic 
sub-programs that make-up the total algorithm. Finally, a 
detailed description of each program will be provided, 
emphasizing the method of implementation of the numerical 
solutions, where applicable, and giving details of a 
procedure to run the individual programs. 

B. THE TEKTRONIX 4052 GRAPHIC COMPUTING SYSTEM 

The host processor used in the TDRL is the Tektronics 
4052 Graphic Computing system. It is a versatile, stand 
alone microprocessor. As such, its capabilities and 
limitations are important driving functions that shape the 
structure of the Operating System. 
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The 4052 has three memories. They are; 

1. Read Only Memory (RDM) — a permanent memory that 
contains the System's intelligence. 

2. Random Access Memory (RAM) --worlcing memory of the 
processor's CP3 . 

3. Line Buffer — small temporary memory that allows data 
to be written and edited on a display screen prior to 
releasing it to the RAM. An important editing 
feature. 

Built-in peripherals for the Graphic System (GS) are a 
smart keyboard (primary input) , a CRT visual display screen 
(primary output) and a magnetic tape unit (mass storage). 
Three extenal peripheral devices extend the versatility of 
the system. They are a Tektronix 4631 Hard Copy Unit (makes 
paper copies of display information) , a 4924 Digital 
Cartridge Tape Drive (additional mass storage) and a Digital 
Processing Oscilloscope (data acquisition). A block diagram 
for the system is shown in Figure 4.1. [Ref. 18] 

The Processor is the main computing device. It 
corresponds to the Central Processing Unit (CPU) found in 
larger systems. It maintains the "firmware" that allows the 
GS to direct system operations, decode instructions and 
perform arithmetic and logic operations. 
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Figure 4.1. Block Diagram of TDRL Graphic System 

The General Purpose Interface provides the means for the 
system to communicate information to or receive data from 
any external peripheral that is compatible with I2EE 
Standard 488-1975. This interface is a bit-parallel, 
byte-serial system capable of transfering ASCII code or 
machine dependant binary code up to 8— bits wide at a maximum 
rate of 250-lc bytes/second. 
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•\ magnetic tape data cartridge is used in the internal 
magnetic and peripheral tape drive units to provide up to 
300-iC bytes of additional mass storage per tape. Programs 
and data may be stored on tape in specifically defined 
"files'*. Piles may range in size from 768-bytes up to the 
limit of the magnetic tape storage capability. Files may be 
either ASCII or binary in format. All files in the 
Operating System are stored in Binary format to take 
advantage of a peripheral "firmware" package called the 
binary loader. Binary data transfers between the processor 
compiler and magnetic tapes are faster than when stored in 
ASCII. This is because the processor uses binary format for 
internal processing and the conversion from ASCII to binary 
format is eliminated when data is stored binary. 
Additionally, a stored binary program can be transfered to 
the processor memory from a specified peripheral device 
without disturbing variables and associated values 
previously defined. This feature is not available to ASCII 
programs . 

The GS permanent ROM has been programmed to respond to 
BASIC (Beginner's All-Purpose Symbolic Instruction Code), a 
high-level programming language first developed at Dartmouth 
College. The Operating System has been correspondingly 
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written in this language. However, BASIC as used by the GS 
has been extended to fit the language to the specialized 
capabilities of Graphic System. These extensions are in the 
areas of graphics, file system access, unified handling of 
input/output operations, matrices, character string 
manipulation, high-level language interrupt handling, and 
operating system facilities. A key difference between 
extended GS BASIC and most orher BASIC languages is that 
most keywords and their parameters can be evaluated 
independently of program control. The result is a rich, 
versatile, yet simple to use programming tool. 

Three peripheral ROM firmware packages have been added 
to the TDRL system to tailor it to specific needs of the 
laboratory. They are the Signal Processing ROM Pack No. 1, 
Signal Processing ROM Pack No. 2 (F?T) , and the EDITOR. A 
description of each is provided in Appendix A. 

C. OVERVIEW OP THE TDRL OPERATING SYSTEM 

The TDRL Operating System has been specifically designed 
to be higly operator-processor interactive. This provides 
the greatest degree of flexibility in the acquisition, 
processing and evaluation of electromagnetic time domain 
transients for target imaging. 
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The total Operating System occupies 83318-byres of mass 
storage on TDRL Library Tape 1 in three files. It requires 
an additional 201628-bytes of storage in 62-files on TDRL 
Library Tape 2 (Data Storage) . The Operating System is 
logically divided into three-individual programs: 

1. INPUT — for the acquisition and averaging of signals 
of interest. 

2. MATH — performs the significant mathmatical processing 
of the acquired signals for inverse scattering. 

3. GRAPH — provides a visual display to the operator of 
the results of the processing. 

TABLE III 

Mass Storage TDRL Programs 



PROG HAM 


FILE # 


BYTES 


% 


INPUT 


2 


24552 


12 


HATH 


3 


14256 


7 


GRAPH 


4 


42624 

1 


22 



Mass storage requirements of the Operating System 
programs are given in TABLE III. The storage requirements 
listed in TABLE III are for the programs only and do not 
include requirements for the necessary definition of arrays, 
strings or numeric constants and variables. 
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All programs are designed to be as nearly independant of 
each other as possible, after an initial run, to allow for 
maximum flexibility and utility in data usage and 
evaluation. The most dependant program is MATH which cannot 
be used without inputs from either INPOT or GRAPH. The most 
independant program is GRAPH which is capable of being run 
with inputs from mass storage only. 

A general flow diagram of the operating system is given 
in Figure 4.2. Specific named subroutines have not been 
shown. The next section will describe in detail the 
individual programs of the Operating System and their 
specific requirements and processing methods, 

D. OPERATING SYSTEM PROGRAM DESCRIPTION 

This section will describe in detail the individual 
programs and their subroutines that make-up the Operating 
System for the TDRL. A general description of the purpose 
of each program will be provided first. Sources of inputs 
and outputs will be noted, A general flow diagram for the 
specific program will be provided sho-wing how all 
subroutines interact. This general overview will then be 
followed by a more detailed description of each subroutine 
in the program. The purpose of the subroutine will be 
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Figure 4.2. General Flow Diagram of TDRL Operating System 



noted. Soarcas of input and output will be defined. A 
discussion of applicable algorithms and the flow of 
information will be also provided. Examples of visual 
displays applicable, potential problems and correction 
methods, will be discussed when warranted. 

As noted, the following commentary will make frequent 
reference to program generated messages. Examples of these 
messages will not be given in the body of this text, but can 
be found in Appendix F, a general run of the Operating 
System. 

1. Program INPaT 

INPOT provides the means to initialize and drive the 
total Operating System. This sub-program is stored on 
Library Tape 1, File 2. Memory space required for its 
implementation in the processor CPU is 24552-bytes. With 
all arrays and constants defined, a total of 42469-bytes of 
memory space may be required. A listing of INPOT can be 
reviewed in Appendix C. 

INPOT is the front-end program for the Operating 
System and as such, will generally be the first program to 
be initiated for TDRL measurements. Program start-up is 
begun by placing Tape 1 into the 4052 built-in magnetic tape 
drive, and placing Data Storage Tape 2 in the 4924 Remote 
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Tape Drive. INPOT is then normally placed into the CPO’s 
working memory under operator control by the issue of the 
following commands from the keyboard: 

FIND 2 
CALL "BOLD" 

Program execution is next begun by typing the command ’’RON” 
from the keyboard. 

If this is the first attempt to run the Operating 
System in conjunction with the Digital Processing 
Oscilloscope (it is assumed the DPO is "OS'* whenever INPOT 
is employed ) , the following message may appear on the GS CRT 
following the issuance of "FIND 2": 

"So SRQ on unit in immediate line — message number 43," 

If the program is being run from any line prior to line 140 
(the usual case), this message may be ignored and the 
program will automatically correct the error message and 
continue processing. If the program is being run subseguent 
to line 140 (the case when returning to ISPOT from some 
other internal file) , the following command must be entered 
from the keyboard: 

POLL N,i!;1 

Following this, the regular sequence of commands may then be 
issued to start the system. 
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As noted, INPOT aay alternately be loaded and 
initiated under program control from GRAPH (File 4) . 
Execution will begin following line 180 in this case. This 
return to INPOT occurs in normal Operating System (OS) 
operation after an initial run has been completed and 
information is then desired for an entirely new target. 

The program specifically provides for the means to 
input information concerning the transient electromagnetic 
response of a target. It acquires, averages and stores all 
pertinent waveforms. INPOT also allows for the easy storage 
of specific antenna and target dimension parameters for 
comparison when the target is known exactly. This data is 
permanently maintained in a library on Tape 2 in F'iles 5-49. 
INPOT further initializes most of the parameters used later 
in Program GRAPH for the visual display of the processed 
results. 

Data inputs are from the keyboard or internal 
memory. Outputs are to internal memory, mass storage and 
Pile 3 (HATH) . INPOT is highly operator-processor 
interactive. Numerous directives and visual displays guide 
the user. 

The structure of the program consists of 
11-subroutines. One (DRIVER) is of general application to 
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the program. Five (ICAVG, INWAV, DIWAV, AUGWAV, SCATWA7) 

are specifically related to input, acquisition, averaging 
and storage of relevant information concerning the 
time- domain transient response of the target in question. 
The remaining five subroutines (ANTLOC, TGTLOC, NSWPAR, 
ANTIil, TGIIN) provide special services for the operating 
system. A general flow diagram for input is provided in 
Figure 4.3 A description of each subroutine follows, 
a. DRI7ER 

This is the main driving routine for INPUT. It 
consists of 129 steps (40 S of the program) from line 100 
through line 1260. It is initiated directly under operator 
control at line 100, or under program control from GRAPH in 
line 180. Data is input from the keyboard or from the mass 
storage files located on Tape 2. DRIVER will pass program 
control directly to MATH when directed. 

Lines 100-300 initialize various variables and 
dimension various arrays. Line 150 clears the General 
Purpose Interface 3us (GPIB) of the DPO service request. 
Flag ?7 (indicates whether the acquired data is destined for 
use in Time Domain or Prony Methods) is set to 0 (Time 
Domain) in line 160. P8, the run number, is initialized to 

zero in line 170 and will be incremented by one in line 210 
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Figure 4.3. General Flow Diagram of INPUT 



for every run (including the first) . 7arious string 
constants are initialized in line 240-250 and 280-290. Flag 
?9 (indicates the need for initializing input parameters) is 
set to 1 (input parameters reguired) . Lines 320-530 provide 
a menu of subroutines and functions which are available for 
Operator selection. The selection of a process from the 
menu is input from the keyboard into system variable PO. An 
example of the menu is given in Appendix F. 

Lines 550-570 direct the program flow to the 
appropriate portion of DRIVER. If the operator has selected 
one of the first four subroutines, and this is an initial 
run of the OS, then the time for a full screen sweep of the 
DPO must be inputted to the program from the keyboard. The 
sweep rate is put in in seconds in T1 in line 640 and must 
be greater than or equal to the time reguired for the the 
DPO beam to sweep completely across the DPO screen. 

Program control is now passed to line 850-1070 
where the following run descriptive parameters are input; 

1. Date 

2. The distance the target is from the transmitting 
antenna 

3. The antenna number selected 

4. The target number selected 
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5. Explanatory remarks 

These are all keyboard entries. An example of the messages 
that query entry of the above parameters is given in 
Appendix F. 

If the run is not the first for the Operating 
System, a different series of inputs are queried. In this 
case, program control is passed to lines 670-840 which 
generate messages on the GS asking the operator if there is 
a change in the original location of the target, if a new 
antenna is to be used, or if a defined target is to be 
changed. The operator inputs the alpha character *'Y” or "N” 
for "yes” or "no" from the keyboard to reflect the current 
situation. Appendix F displays these messages. 

Lines 1000-1040 generate the input of the 
antenna parameters. These parameters are drawn from files 
stored on Data Tape 2. Each antenna file, located in one of 
the files numbered 5-24, is a 63-element array, returned in 
A5. In elements 1 through 60 are located: 

1. the run number for this particular antenna- target 
combination 

2. the antenna length compensating parameter for the 
antenna-target combination 



108 



3. the antenna amplitude compensator for the 
ant enna- target combination 

Is can be seen, antenna parameters are actually three 
elements within the larger array, each three elements 
identifying a unigue antenna-target combination. The 
parameters in elements 51-63 are system antenna parameters 
used for unknown targets. They are simply the sum of all 
run, antenna length compensator updates and antenna 
amplitude compensator updates. The average of the last two 
elements is applied to signals acquired if the target is not 
defined from the target library. This is shown in lines 
1024-1025. 

The antenna amplitude compensator is a very 
sensitive function of the distance the target is from the 
transmitter due to EM propagation attenuation. In line 
1025, for unknown targets, and in line 4152 for defined 
targets, these losses are taken in consideration by finding 
the square of the ratio of the targets true location on the 
plana with reference to the transmitting antenna and the 
control distance (1.27m). 

After the definition of variables is complete, 
DRIVER passes control to the appropriate function in line 
1090. Line 1100 returns the display to the Menu for further 
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function selection or to line 1170 where the DRIVER passes 
control to MATH and terminates any further INPUT processing. 

In summary, INPUT DRIVER reguires input of 
selected program routines and initial parameters by the 
operator from the keyboard, and requires the OS to determine 
whether the current run is the initial for the day or a 
subsequent one. The parameters are saved for later graph 
requiremts. Funtion selection drives the appropriate 
subroutines and routines. A flow diagram for GRAPH DRIVER 
is provided in Figure 4.4. 

b. Subroutine AHTLOC/TGTLOC 

ANTLOC and TGTLOC are similiar subroutines in 
that they locate the appropriate antenna and target storage 
arrays filed on Tape 2. Antenna parameters are located in 
Piles 5-24. Target parameters are located in files 26-45. 
Lines 4000-4060 contain ANTLOC. TGTLOC is found in lines 
4100-4160. TGTLOC has the additional function of defining 
the system antenna parameters for the current run, as shown 
in lines 4149-4153. Thera are no operator-procecessor 
’interactions during execution of these subroutines. 

c. Subroutine DIWAv/INHAV/AOGWAV 

DIHAV, INWAV and AU3HAV are similar subroutines. 
They drive ACAVG to ioput the direct waveform, the incident 
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waveform, and the augmented waveform. They allow the 
selection by the Operator of the number of signals to be 
acquired for each waveform by the DPO (total composite 
waveforms from the sampler) and averaged by the OS. Bach of 
these subroutines is accompanied by a message on CRT 
describing the procedure to be employed in acquiring the 
particular signal. These messages are reproduced in 
Appendix F. DIWAV is located at lines 3000-3170. It saves 
data in File 50 of Tape 2. INWAV is located at lines 
3200-3380. It saves data in File 51. AOGWAV is located at 
lines 3400-3590. It saves data in File 52. 
d. Subroutine ACAVG 

ACA7G acquires and averages the appropriate 
waveform interactively requested by the operator from the 
processor keyboard. It also removes any dc signal that 
might be included in the acquisition. This feature allows 
the operator to position the desired signal anywhere on the 
DPO screen and still maintain a true zero axis in later 
graphic displays. Note, the signal waveform must be 

completely within the DPO screen area. ACAVG is 74-steps 

long { 23 % of the program) and is located in statements 
2000-2730. It is driven by one of the appropriate input 
waveform subroutines previously described. It returns 
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results to that subroutine for later storage. It also 
displays system information for the current signal being 
acquired to the DPO screen. 

Lines 2990-2130 initialize ACAVG by requesting 
the operator to input the number of waveform signals to be 
averaged in N2. This value is then returned for 

verification. In lines 2140-2540, the OS writes the current 
acquisiton number of the waveform to the DPO screen, 
acquires and stores the signal into DPO memory and then 
transfers the signal to the 4052 for further processing. 

Lines 2170-2370 perform the count function. 
Depending on the number of waveforms desired, the count will 
be for every signal (N2 lass than 21) , for every five 
signals (N2 less than 101) or for every ten signals (N2 
greater than or equal to 101). This is done to reduce the 
time required for acquisition of all signals. 

Lines 2410-2430 instruct the DPO to store and 
hold the waveform in DPO memory. Line 2420 stops processing 
to allow for the DPO to fully acquire the waveform in a 
single sweep. T1 is the length of the wait. 

In lines 2440-2450, the GS instructs the DPO to 
transfer the digitized signal from DPO memory to the 
microprocessor temporary storage array XO. Lines 2470-2500 
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then search the first eight values of the array for any 
obviously out of range values input by poor DPO blanking on 
the signal front end. If such a value is found, the array 
element is set to discrete value 511, corresponding to a 
scaled value of 0, The waveform is next added to the sum of 
the previous waveforms for the current acquisition in line 
2530. The resultant matrix is the time sampled values of 
the signal in BO: 

N “■ 

BO = Z XO. (4.1) 
i=l ^ 

If further acguistions are necessary, the entire process 
described above is repeated. If no further acquisitons are 
to be done on the current signal, the ensemble average time 
sampling matrix XO is determined in line 2550: 



XO = B0/N2 



(4.2) 



4ny dc offset is eliminated in line 2560. 

It is interesting to note the time required to 
acquire and average signals using the solid state OHF pulse 
generator. For ten averages of a signal, 42 seconds of 
elapsed time is needed, or about 4.2 seconds per acquistion. 
Of this, about 1 second is needed to write information to 
the DPO screen per acquisition. For 21 averages, about 65 
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is required, or about 3,1 



seconds of acquisition time 
seconds per acquisition. This is about a 35 % reduction in 
acquisition time overall per signal. The improved time per 
acquisition is due to the reduced frequency with which data 
in the form of the wavecount is written to the DPO screen. 
For a relatively small increased investment of time, a 
tremendous advantage is realized in averaging signals and 
reducing polluting noise. A further decrease in acquisition 
time per signal is realized whan 101 or more averages are 
specified. However, as will be discussed in Chapter 7, lack 
of total signal stability may produce errorneous results for 
such a large number of acquisitions due to the length of 
time required. 

In lines 2570-2820, the waveform scale factors 
are determined and transferred to the processor. The values 
are returned in (vertical magnitude scale factor) and 
(horizontal time scale factor) . The discrete, digitized 
verticle values are appropriately scaled to the correct real 
time values in line 2700. The time between a single 
discrete sample on the horizontal axis (sampling rate) is 
returned in Z1 at line 2820. 

In summary, ACAVG does the physical acquisition, 
averaging and scaling of each signal. It returns the 
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current waveform count to the DPO CRT for display, and the 
appropriately scaled waveform in XO to the driving 
subroutine for storage in the appropriate file. It further 
returns the appropriate horizontal and vertical axis labling 
scale factors in T and M respectively. The flow diagram 
for ACAVG is provided in Figure 4.5 



The 


basis 


for 


this subroutine 


was derived 


by 


Capt. Hammond. 


[Ref. 


19] 


His program in 


its entirety 


has 



been extensively modified to fit the current requirements of 
the TDHL OS. 

e. Subroutine NEWPAR/ANTIN/TGTIN 

NEWPAR, ANTIN and TGTIN are interactive 

subroutines which allow an operator to input new parameters 
for defined antennas and targets into the mass storage 
library on Data Storage Tape 2. 

NEHPAR is the driving routine for AHTIN and 
TGTIN. It receives its input from the DRIVER and selects 
either ANTIN or TGTIN. It returns the new parameters to the 
appropriate file on the storage tape. NEWPAR then determines 
if further inputs are required. If not, program control is 
returned to DRIVER. 

ANTIN defines the values of the parameters of 
the new antenna. It requests the operator to input the 
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antenna number. An array of 63 elements is then generated, 
each third element being filled with an initial antenna run 
number equal to 1. All other elements are sat to zero. If 
the operator specifies an antenna value greater than 20, the 
following message is returned. 

••STORAGE AREA LIHITED TO 20 ASTEHHAS. WAHT TO COHTIHUE?" 
This is necessary as current space limitations limit the 
number of antennas to be defined. 

TGTIN defines the values of the parameters of 
new targets. It requests the operator to input the new 
target number as well as the length and radius of the 
target. The carraige return must be pressed after each 
entry. 

NEHPAS, ANTIS and TGTIN are located in lines 
4200-4800. They comprise ^^% of the program, 
f. Subroutine SCATWA7 

SCATWAV determines the target backscattered 
waveform based on the difference of the augmented and 
incident waveforms. It returns the backscattered waveform 
to storage at file 53 on Tape 2. It also sets flag P7 to the 
appropriate value required for either Prony processing 
(P7 = 1) of Time Domain processing (p7=0) . Program control 
is then returned to INPOT DRIVER. 
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2 * P r qg raia .'lATH 

MATH parforas aost of tha significant mathaatical 
computations for tha OS. Its correct performance is central 
in the production of accurately processed signals. MATH is 
stored on Tape 1, File 3. It requires V4544-bytes of CPH 
memory space. With all arrays and constants defined, a 
maximum total of 49108-bytes of memory are required. 

This sub-program is initiated under program control 
only from INPUT and GRAPH, There is no interaction between 
operator and system during MATH'S execution. 

Inputs to MATH are from INPUT, GRAPH or Tape 2. 
Values are returned to mass strorage or GRAPH. 

MATH performs the following functions for the OS: 

1, Fast Fourier transforms the Direct and Backscaztered 
waveforms . 

2. Determines the frequency domain impulse response, 

3, Performs inverse Fourier transform to find the time 
domain inpulse response, 

4. Determines the step and ramp response for the target. 

5. Determines the physical optics shape of the target, 

6, Optimally filters the frequency domain impulse 

response, yielding a smoothed time domain impulse 
response, ramp response and physical optics shape of 
the target under noisy conditions. 
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7. Compensates the ramp response and optics shape for EM 
propagation attenuation. 

8, Updates optimal antenna parameters for the antenna 
being used, based on the targets ramp response, 
automatically. 

HATH is composed of four subroutines and DRIVER. A 
description of each follows. A general flow diagram for HATH 
is gi7en in Figure 4.5. 

a. MATH DRIVER 

The driver routine for MATH is located between 
lines 100 to 1210 and comprises 59% of the total program. 
In addition to initializing the program and driving the 
subroutines, DRIVER performs most of the significant 

mathematical computations in straight line data flow 

fashion. 

MATH DRIVER is initiated at line 100 from INPUT 
or line 130 from GRAPH. Parametric data is input directly 
from INPUT. It is recovered from mass storage in line 160 
if GRAPH initiates MATH. Lines 100-221 serve to initialize 
the local program parameters. 

Line 250 inputs into XO from mass storage file 
53 the signal (backsca tter ed waveform) which represents the 
convolved transient response of the transmitting antenna. 
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the receiving antenna, and the target. This data is then 
used to drive subroutine FFT. XO is returned from FFT in M 
and P as the frequency domain magnitude and phase and saved 
in arrays E and S3 respectively. 

A similar procedure is performed in lines 
290-340 for the convolved transient response of the 
transmitting and receiving antennas without target on the 
ground plane. The magnitude and phase are returned in 
arrays F and P3 respectively. 

In lines 360-420, the frequency domain impulse 
response of the target and noise optimization are performed. 
Magnitude of the impulse response is returned in EO, with 
phase returned in FO. The algorithm for determining the 
impulse response and performing optimization were developed 
in chapter III and will not be repeated hera. 

Lines 440-540 perform the inverse Fourier 
transform of the frequency domain impulse response of the 
target to find its time domain equivalent. The frequency 
domain data entering at line 440 is in polar form. This 
form is most suitable for the previously described 
mathematical manipulations. It is not a suitable form for 
the inverse Fourier transform method employed by the 4052 
system. In order to perform this operation, the magnitude 
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and phase arrays must first be converted into two source 
arrays containing real and imaginary data. Next, the first 
element of the real array must contain the signal dc term, 
while the last real element must be the Nyguist frequency 
value. As both these elements represent pure real numbers, 
the first and last elements of the imaginary source array 
muse be correspondingly zero. 

Lines 440-470 establish a FOR/HEX? loop that 
returns real terms to M and imaginary terms to P. Real 
terms are established by multiplying the polar magnitude, 
EO , by the cosine of the corresponding phase FO: 

M(j) = E0(j) X (cos (F0(j)) ( 4 . 3 ) 

Imaginary terms are established by multiplying the polar 
magnitude by the sine of the corresponding phase term: 

P(j) = E0(j) X sin(F0(j)) C^.4) 

INPUT had previously removed all dc levels prior 
to NATH, so N(1) and 3(257), the dc and Nyquist terms, are 
set equal to zero in lines 480-490. Lines 500-510 set the 
corresponding imaginary terms to zero. Note that both of 
these arrays must be 257 data points long, as the output 
array will be 512 data points. 
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Data 



is now correctly formatted in the two 



source arrays so that the command "INLZAV” can interleave 
the real and imaginary data into a single destination array, 
lO , whose format is acceptable to the GS command. 

Line 540 calls the system inverse Fourier 
transform command which transforms the spectral data 
into waveform (time- doma in> data. The spectral data 
originally in 10 is overwritten by the new time domain data 
and returned in 10. 10 is now the time domain impulse 
response. 

The inverse Fourier transform can be expressed 
mathematically by the following summation: 



X(n) 



(1/N) 



N-1 

Z 

k=0 



X^(k) 



for n=0,l, . . . ,N-1 (4.5) 



In the above equation, N refers to length of the array 
argument, and n is an index used in generating the various 
elements of the output array. X (k) is the k-th complex 
Fourier coefficient, and X (n) is the n + 1-th element of the 
real data output. [ Eef . 20] The 4052 uses an extremely 

fast method for IFT computations, the Sande-Tukey 

decimalization- in- frequency algorithm , 
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Line 620-720 save the time domain impulse 

response (file 54) , and frequency domain direct waveform 

(file 55) , the backs catt ered waveform (file 56) and the 
impulse response (file 57) in mass storage. 

Next, the time domain step response of the 

impulse is returned to XO in line 740 by using the 4052 

routine. This is an integration calculated using the 
trapezoidal rule for approximating the definite integral as 
follows : 

B = 0 

B (t) = B(t-1) + .5 (A(t-1) + A(t)) Cm 6^ 

for t = 2, 3 , . . . ,N 

where: A = the source array 

B = the destination array (result) 

N = the number of elements in the array 
Any dc value in the result is removed through the actions of 
lines 750-76 0. 

Next, the time domain ramp response is 
determined by integrating the step response in XO. The ramp 
is returned in HO. The HATH DRIVES outputs the time domain 
ramp and its minimum value in R2 to Subroutine OPTICS which 
returns the physical optics shape (shape) and ramp response 
compensated for electromagnetic propagation attenuation. 
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Note the shape does not as yet have the correct amplitude 
values in its array elements. 

If the run is the first for the target (prior to 
optimization) and the target is one defined in mass storage, 
MATH DBI7ER inputs the last computed antenna parameters for 
the compensator for shape length (S5) and radius (S6 ) , as 
well as the actual defined langhth (C6) and radius (C5) to 
subroutine ANTDP. S5 and S6 are then automatically updated. 
The physical optics array is then multiplied by the 
compensator value S5 so that its elements reflect the 
estimated dimensions of the targets in inches. This is 
based on the experimentally observed fact that a particular 
size and configuration of stationary target will return some 
baclcscattered field relative to its size and shape. That 
this value is not exactly the same from one run to the next, 
but only approximately so, is due to certain stastically 
determinable effects produced by noise, small differences in 
orientation, etc. A large number of runs should produce an 
averaging value which takes into account these factors. If 
the target is unknown (the general case) , or the current run 
is for noise reduction, previously computed values of the 
antenna parameters are recalled from mass storage. For the 
unknown target these compensators are the averages for all 
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targets for all runs for a particular antenna. In fact, 
OPTICS is not part of the data processing flow, the 
compensators having been previously determined in INPUT at 
lines 1024-1025. 

Following a known target run, the minimum value 
of BO is returned to B8. If the run is not for optimization 
(L1=1) , then B6, an amplitude compensation optimization 
factor is set equal to B8 and the actual length of the shape 
is retained in Q9. For graphing purposes, it is desirable 
that B6 and Q9 not vary from optimization run to 
optimization run. MATH DRIVER will input BO, Q9, and B6 to 
subroutine OPTIM, which compensates the amplitude and length 
of the physical optics shape for the low pass filtering done 
to reduce noise. An arbitrary target is treated in exactly 
the same manner. (Note, this effects the shape ONLY and not 
the time domain ramp response during optimization.) 

The above system has worked remarkably well in 
deducing the actual physical parameters of both known and 
unknown targets. For known targets, accuracy is invariably 
^ 0 % or less. For unknown targets, this accuracy 
deteriorates to about 15% or less. If the data is treated 
in the raw form, with no compensation applied, accuracy is 
reduced to 25% or worse, up to 50% for the unknown target. 
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This can be noted on the graphic output of the shape where 
the compensated, uncompensated and known (when given) values 
of length ace compared run to run. 

Lines 1000-1140 save the updated antenna 
parameters in the appropriate files, the shape array in file 
58, the time domain ramp response in file 59, and 
significant system parameters in file 60, all on mass data 
storage Tape 2. 

Lines 1160-1190 is MATH DRIVES ending routine, 
deleting variables and passing control of the program to 
GRAPH (file 4) on Tape 1. 

In summary, MATH DRIVER is the principal vehicle 
for OS processing for inverse scattering. A flow diagram is 
given in Figure 4.7. 

b. Subroutine FFT 

FFT is the fast Fourier transform subroutine 
used tc convert the time domain Direct and Backscattered 
waveforms to spectral data. This subroutine is located in 
lines 2000-2140 and requires 7S of the total program space. 

XO contains the real-valued signal data on which 
the Fourier transform is to be performed. It is input from 
MATH DRIVER. Line 2060 calls the GS command '«ffT" and 
performs a fast calculation (about 4 seconds) of the DFT 
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Fissure 4,7. Flow Diasjrapi--MATH DRIVER 



(discrete Fourier Transform) , 
the summation: [Ref. 21] 



expressed mathematically by 



Xd<k) 



n=0 



for k=0,l, , . . ,N/2 



(4.7) 



where: N = array length 

k = index used in generating various 
Fourier coefficients 

Xj(k) = Ic-th Fourier coefficient 

X(n) = n + 1-th elemant of the real data 
input 

After the command has executed, the original 
real data will be overwritten by the results of the FFT 
computation. The output data is in complex rectangular fora 
with the first element of the array XO containing a real 
number representing the dc value of the signal, the second 
element containing a real number representing the value of 
the discrete Fourier transform at the Nyguist frequency, and 
the remaining array elements containing alternately real and 
imaginary Fourier coefficients. 

This arrangement is not suitable for the 
division that later takes place. In line 2070, the 
rectangular array XO is converted to polar form by calling 
the GS command ’’POLAR”. Magnitude components are returned 
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in array H and phase components are returned in array P. At 
this point, subroutine FFT returns control to MATH DRIVES 
with the arrays XO, M and P for further processing, 
c. Subroutine OPTIC 

Subroutine OPTIC manipulates the ramp response 
of the target to recover the physical optics shape, but it 
does so without effecting in anyway the raw data. It later 
compensates both the ramp response and the shape for EM 
propagation attenuation experienced along the body of the 
target. Subroutine OPTIC resides in lines 2200-2700 and 
comprises 25% of Program MATH. 

The ramp response RO and the location of its 
minimum value within the array are input into OPTIC from 
MATH DRIVES. The mean value of the response is calculated 
in lines 2250-2260 and returned in R3. This value is used 
later in determining the endpoints of the shape array. 

Lines 2290-2500 is a FOR/NEXT loop that 

seperates out the targets shape from the ramp array data. 

The search for valid data begins at the point in the array 
equal to the location cf the minimum value. It is assumed 
that the target will always return a response larger than 
any that might be due to range noise or other excitation. 

For the size of the targets being used, and the the 
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configuration of the present source, this has always proven 
to be an accurate assumption. The FOR/NSXT loop searchs 
first for a trailing-edge endpoint. This is done in 
ascending time steps in search of the point where the array 
elements become zero or positive, representing a zero-axis 
crossing. For every value of SO less than zero, a 
corresponding value is input in BO, the shape array. If 
30 (k) is a positive value, representing a zero-crossing, the 
F03/NEX7 loop is terminated. If BO (k) is negative, OPTIC 
then checks the location of the point. It is important to 
do this as noise and instability of input signals from an 
acquisition- to-acquisition run may erroneously shift-zero 
crossings or preclude zero-crossings altogether. In these 
cases, determination of the shape must be somewhat modified 
as described below. 

In line 2350, two situations are checked. 
First, is the preceding array element of RO more negative 
than the present element and is the present element greater 
that 90% of the ramp response average value. If 
affirmative, the next array element is fetched. A positive 
result of this check indicates that datum is on the skirt of 
the ramp response, and that the data magnitudes are 
progressing towards a zero-axis crossing normally. If the 
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checJc returns as negative, that is, the present array value 
is nore negative than the preceding element, than an anomaly 
has occurred either due to noise, or time-axis shifting. In 
this case a special endpoint routine must be employed. 

The second check is to determine if the array 
element is on the main signal portion of the response and 
not the skirts. If the element amplitude is greater that 
90% of the mean value, than it is assumed that the element 
is on on the skirt ir regardless of the relationship between 
the amplitude of the present and preceding elements. The 
next datum is fetched. It can be seen that this arrangement 
is subject to some error if a noisy element is encountered 
just at the point where the ramp response passes the 90-th 
percentile of the signals mean value. This error is 
considered relatively minor. 

The 90-th percentile was chosen following 
numerous observations that noise or time-shift in the main 
ramp response normally do not produce fluctuations that 
cross this boundary. Zrror in the shape representation will 
be introduced if it does, in fact, cross this level. The 
result will be a fore- shortened drawing of the target. 

If the need for an endpoint routine is 
indicated, processing goes to lines 2360-2U70. In line 
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2360“2380, the distance the last valid data element 
(30(;c-1)) is from the zero axis is determined. In line 
2360, the distance between the previous two valid data 
elements is found and returned in W (3) : 

W{3) = R0(k-2) - RO(k-l) (m . **) 

The above is for the ascending time series run of the 

FOR/NEXT loop (trailing edge) . For the descending time 
series, the values in the parenthesis have opposite sign. 
This difference represents an approximation of the slope of 
the skirt. A problem with this approximation is that it 
often occurs where the slope is not as great as the average. 
The result is a somewhat elongated figure. It may be 

possible to defeat this problem by taking more than two 
points for the average. 

In lines 2370-2380, a check is made to see if 
the last valid data point is within 10-times or less of the 

difference determined above. If it is, the actual 

difference is returned in Q1 and used to drive the values 
placed in BO to a zero-crossing. If the difference is less 
that 10% of the amplitude needed to effect a zero-crossing, 
then one- tenth of the actual distance is returned in Q1. In 
this manner, the shape rapidly goes to zero and an endpoint 
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is detenainad close to the estimated real endpoint had a 
natural zero-crossing occured. The endpoint is located with 
a Foa/NEXT loop between lines 2U20-2470. This FOa/NEXT loop 
simply iterates the values of BO by Q1 until a zero-crossing 
is effected. Whenever BO becomes positive, all loop 
processing is terminated and the endpoint is returned in 
W(1) for the trailing skirt and W (2) for the leading skirt. 

If the above process was for the leading edge of 
the ramp response, the loop is reentered and the trailing 
edge check of the data is performed. The method employed is 
identical with the above process for the leading edge. Note 
that up to this point, the ramp response in RO has been 
totally unaffected by any processing in OPTICS. It has been 
used simply as source of data. 

Once the target response has been located in RO 
and entered in BO, then both arrays can be compensated for 
propagation losses experienced. This is done in lines 
2570-2640. The value of the nominal distance the target is 
from the transmitting antenna is returned in Z6. A FOR/NEXT 
loop in lines 2590-2640 corrects the arrays for the 
estimated propagation losses. Line 2500 computes the 
distance along the body the wave propagates for each sample 
point in RO: 
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H = (Q - W(1) * Z1 * (10*8) n 



Cm -8) 

where; (Q - W(1)) * Z1 = twice the time of travel 

of the propagated wave 

In line 2610, the propagation loss compensation factor is 
returned using; 

u = (4.10) 

2Z^ 

where; Z = range from transmitting antenna to 
target face 

It is also assumed that the distance between transmitting 
and receiving antenna is Z. Z is in fact the measured value 
of distance read from the "TIME/DISTANCE” dial on the DPO. 
Each value of RO and BO is then multiplied by this 
compensation factor in turn. Processing continues once all 
values for the computed length have been appropriately 
amended. 

Line 2260 is a FOS/IIEXT loop that transfers the 
values in the array BO to the beginning of the array to 
facilitate graphing. Once this is accomplished, program 
control is returned to MATH DRIVER with the shape unsealed 
for actual dimensions of the target (values are equal in 
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amplitude to the compensated voltage magnitudes of the ramp 
response) . The ramp resonss is returned uneffected by any 
of the above processing except for propagation attenuation 
compensations. Figure 4.8 provides a flow diagram for 
Subroutine OPTICS. 

d. Subroutine AIITOP 

This surboutine updates the specific antenna 
length and radius compensation values if the target is one 
defined from storage on Tape 2. If the target is arbitrary 
in the general case, no updating occurs. ANTOP is located 
between statements 3000-3080 and comprises 4% of the 
program. 

The unsealed shape array BO is input to ANTOP 
from aATH DRIVER. In line 3060, the amplitude compensating 
and scaling factor is computed using; 

3fe= (A ♦ S*B)/(C*B) 

where: S6 = scaling/compensating amplitude factor 

C = radius of target 

B = minimum magnitude of physical optics 
array 

A = the number of previous times this 
antenna/target combination was used 

S = sum of previous compensation values 
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The result is simply the averaged value of 
compensation/scaling for all runs using this particular 
antenna and target combination. 

Line 3070 performs a similar computation for S5, 
updating the compensation/scaling factor for the target 

length. The formula is identical except that the defined 

target length replaces the target radius and the observed or 
measured length from the plane replaces the minimum 

magnitude of the shape array. The concept is that slight 

differences in target returns due to variations in target 
location and orientation, noise, etc., will be averaged out 
over the long term for any particular antenna. When the 
target is arbitrary and there is no a priori information 
then AWTOP is bypassed and the average of all the 
compensation/scaling factors for all targets for the 

particular antenna is used. This requires that any new 

antenna be ’’seeded'* with compensation/scaling factors from 
at least one run using known targets prior to attempting to 
run an arbitrary target on the antenna. What results is a 
’’poor man’s” compensator for the currents on the body of the 
target based on empiracal observations of many target 
responses. It is not as accurate as the full solution of 
the time domain integral equations would be, nor is the 
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exact target shape returned. What is obtained is an 
excellent estimate of the target, whether known of unknown, 
based on an historical compilation of data for the TDRL 
system. This concept can be generalized to any system using 
components which are nominally fixed, such as large scale 
bistatic radars. The shape itself, at least for simple 
axisymmetric targets, can probably be returned with the 
development of the proper software routines. Perhaps the 
rate at which the ramp response slopes vary may be analyzed 
for this purpose. What is obtained is a method which 
employes a small scale processor requiring relatively simple 
processing to obtain accurate time domain estimates based on 
the ramp response of the target almost in real time, and 
capable of being generalized to almost any system employing 
the proper source. 

The compensation/scaling factors are returned to 
MATH DHI7ER along with program control, BO is unaffected as 
yet by S5 or S6. 

e. Subroutine OPTIM 

OPTIM compensates the shape array for the 
effects of low-pass filtering of the target responses using 
Riad's Optimization technique. Optimization is done to 
remove high-frequency noise terms and, for the shape, to 
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improve the overall shape of the drawing. In reality, 
without compensation, the drawing is less representative of 
the true targets shape when optimized. Low pass filtering 
will lengthen the overall shape and reduce the radius, the 
amount being relative to the value of lambda, the 
optimization variable. OPTIM allows for this and introduces 
factors which will pin the optimized shape to the original 
form while allowing the benefits (reduced noise ripple on 
the body) of low-pass filter noise reduction to be shown. 
OPTIh is located in lines 3100-3180 and represents 4^ of the 
program. 

30, B6, (original minimum of the unoptimized 
shape) , 38 (the present minimum of the shape) and Q9 (the 
original unoptimized target length) are input from MATH 
DRIVER. The amplitude OPTIM compensation is found in line 
3140 by dividing B6 by B8. It is returned in B7. The new 
compensated target length factor is returned in Q7 and is 
found by dividing the present value of the length by the 
original value. Finally the whole physical optics shape 
array is scaled by the new amplitude factor. The new values 
are returned to MATH DRIVER, along with the original factors 
in 36 and Q9 . In fact, after the initial unoptimized run, 
36 and Q9 are constants until a new target is acquired. 



^ • Prog ra m GRAPH 



GRAPH provides a central means to visually display 
the results of processing performed in HATH . The correct 
evaluation and analysis of the acquired processed data is a 
direct function of the effective reproduction of the 
resultant arrays by GRAPH into a form most useful and 
understandable to the operator. GRAPH is located on Tape 1, 
Pile 4, and requires 44712-bytes of memory space. 

Program GRAPH is normally initiated under program 
control from Program HATH, in line 100. It may be initiated 
directly by the operator if the results of a previous run 
have been retained in mass storage. The commands issued by 
the operator are: 

FIND 4 
CALL BOLD 
BON 

Outputs are to the CRT or a return of control to one of the 
other programs that make-up the TDRL OS. No results are 
returned to mass storage. 

GRAPH subroutines may be divided into three 
categories : 

1. DRIVER-- provides for initialization of the program 
and drives subroutines. It also passes control to 
the next program in the OS. 
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2. AXIS and Labling Sabr3Utin6s--establishes the axis of 
all graphs and lables as appropriate. Provides for 
formatting of descriptive parameters, as appropriate 
builds grids for graphics. Six subroutines comprise 
this group. They are AXIS, TICKS, AXLAB, TITLE, 
GRPAR, GRID. 

3. Inverse Scatter ing--Subroutines in this section graph 
the results of processing in INPUT and MATH. They 
reproduce the four basic input waveforms, the Fourier 
transforms of the Direct, Backscatter and Impulse 
Response, the Physical Optics Shape of the target, 
and the time domain ramp response. Four subroutines 
comprise this group. They are GRIN, GRFOUR, GROP, 
and GRAMP. 

A general flow diagram for GRAPH is provided in 
Figure 4.9. A disscussion of all Subroutines in GRAPH 
follows. 

a. Graph DRIVER 

As for INPUT and MATH DRIVER, GRAPH DRIVER is 
the principal initialiser and controiing routine for this 
program, coordinaring the flow of processing. Unlike MATH 
DRIVER, GRAPH DRIVER does not of itself perform any data 
processing. That function is retained for the subroutines 
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GRAPH DRIVER is located in lines 100-1190 and 
comprises 19^ of the total program. GRAPH is completely 
independant of data input from the other programs in the OS. 
All inputs are from mass storage Tape 2 or the GS keyboard. 
All results are returned to CRT display of the 4052. This 
routine is highly operator-processor interactive. Appendix 
? describes a typical run with messages as produced. 

GRAPH DRIVER is always initiated in line 100 
when being entered external to the GRAPH Program. It is 
initiated from line 290 by returns from internal 
subroutines. Lines 100-300 initialize the program and input 
parameter data from mass storage. Lines 320-750 develope 
the correct Menu, with the function/subroutine selected 
input to PO in line 750. 

The Menu that is used in GRAPH lis shown in 
Apendix F. Host of the fiunctions on the menu are 
self-explanatory. A few remarks concerning items 10, 11, 
12, and 13 should suffice to remove any ambiguity. 

Item 10, '’COSTINHE OPTIMIZATION”, initiates the 
Riad low-pass filtering optimization of the impulse response 
in lines 950-1030. Program control with the operator 



designated optimization factor is input to Program HATH. 



Item 11, ’’NEXT TARSET", passes program control 
to INPOT for the acguisition of an entirely new target. 
Parameters are retained and passed along with program 
control . 

Item 12, ’’FINISKED GRAPHING", terminates 
processing. The program can be reinitiated by typing "RON". 

Lines 820-900 select the appropriate subroutine 
and functions. Line 940 returns to the menu for further 
selection for all items except 10,11, 12, 13. 

The next series of subroutines to be discussed 
perform the graphics axis formatting, Tabling and titling. 
These subroutines set the environment that makes the 
information of use to an observer. AXIS, TICKS, AXLAB, and 
TITLE are subroutines based on programs previously developed 
by Tektronix, Inc. [Ref. 22] All have been modified to a 
greater or lesser degree dependant on the needs of the TDRL. 
b. Subroutines AX IS/TICKS/AXLAB 

Subroutine AXIS drives Subroutines TICKS and 
AXLAB. They perform the general axis definition function 
for the output graphics. Subroutine AXIS serves as the 
imput portal for data from the function subroutine that was 
chosen by the operator. The results returned to the CRT of 
the GS are "neat" ticks, Tabled and scaled to the proper 
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values of the input function. Subroutine AXIS occupies 
statements 2000-2250, TICKS is in lines 2300-2610 and AXLAB 
is in statements 2700-3190. These three subroutines occupy 
18% of program space. 

Data input from the function to be graphed are 
the minimum and maximum horizontal values of the viewport in 
VI and V2, and the minimum and maximum vertical values of 
the viewport in V3 and V4. These constants are in Graphic 
Display Onits (GDO) . GDD's are a measure of the resolution 
capability of the 3S screen. The four arguments defined 
above refer to the actual limits for the graphic information 
of the display. The GS will draw no lines outside these 
limits . 

Also input to Subroutine AXIS from the driving 
function subroutine are the minimum and maximum horizontal 
values of the function itself, '^1 and X2, and the minimum 
and maximum vertical values, H3 and M4. These values are in 
User Data Onits (UDO’s). That is, in units that the user 
defines, such as millivolts, seconds, meters, etc. 

Finally, AXIS receives the value of the flag SO. 
This specifies whether the vertical axis (S0=0) or 
horizontal axis (S0=1) is being labled. The magnitude of the 
horizontal scaling function for the input wave, SI (received 
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froa INPUT) , and the definitions in UDO’s for the vertical 
axis and horizontal axis are required. 

The performance of AXIS, TICKS, and AXLAB is 
adequately described in the Tektronix source references. 
[Ref. 23] This report will describe only those statements 
which have been changed or modified to produce the desired 
output for TDRL purposes. 

Principally, the changes effect the manner in 
which the axis is lable with UDO’s in Subroutine AXLAB. 
These changes lie in statements 2860-3190. Three conditions 
must be checked: 

1. Horizontal (S0=0) or vertical axis (S0=1) lables? 

2. Graphics of a Fourier transform? (4<P0<8) 

3. Check amplitude of maximum vertical value (P6) in 
ODO’s. 

TABLE IV defines the relationship between 
statements in Subroutine AXLAB and the functions to be 
graphed. 

c. Subroutine TITLE 

This subroutine positions the title of the 
graphic above the plotted values. It is located in lines 
3200-3310 and comprises 2 % of rhe program. The value of the 
maximum horizontal coordinate M2 and the graphic computed 
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TABLE IV 



Graphic Labling Relationships — GRAPH AXLAB 



GRAPH 


STATEMENTS 


REMARKS 


GRIN 

GRAHP 


HORIZONTAL 

AXIS 


HORIZONTAL 

AXIS 




2880-2890 


S1*10>1000. Lable Axis 




2910 


S1*10<1000. Lable Axis 




30 90 


Position Dimension Lable 




31 10 


Dimension Lable 




31 40-3160 


Print Dimension Lable 




3180-3190 


Return 




VERTICAL 

AXIS 


VERTICAL 

AXIS 




2860 


Start 




2940-2970 


P(6)<1. Lable axis 




2980-2990 


1<P(6)<10. Lable axis 




3000-3010 


P(6)>10. Lable axis 




31 30-3150 


Print Dimension Lable 




3180-3190 


Return 








GRFOOR 


HORIZONTAL 

AXIS 


HORIZONTAL 

AXIS 




2861-2862 






2863-2864 


"Magnitude" Axis Label 




2865-2866 


"Phase" Axis Lable 




3171 


Print Dimension Lable 




3180-3190 


Return 




VERTICLE 

AXIS 


VERTICLE 

AXIS 




2860-3170 


Same as for Horizontal Axis 


1 


3173-3190 

j 


Return 

1 



NOTE: All graphic routines use lines 3020-3040 
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vertical maxiaura P{6) are input by the graphic subroutine. 
The graph title centered over the graph is returned to the 
CRT. 

d. Subroutine GRPAE 

GRPAR graphs the parameter input by INPDT and 
developed by MATH. This subroutine is located in lines 
3400-3870 and comprises 8S of the program. All graphs have 
the run number, date, target distance from transmitter, 
antenna number, target number, averages taken, optimization 
and any specific remarks printed. The specialized 
parameters printed for particular graphs alone with the 
statement number applicable are given in TABLE V. 

TABLE 7 

Specialized Parameter Prints — GRAPH GRPAR 



GRAPHIC 


PARAMETERS 


LINES 


GRIN 

GRAMP 




3510-3600 


MAXIMOH PEAK VALUE 


aiN-iaOM PEAK VALUE 


RMS VALUE 


MEAN VALUE 




GROP 

1 




3710-3750 


DIAMETER (INCHES) 

- - - - 


;?IDTH (INCHES) 
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Subroutine GRID 



This simple subroutine draws referencing grid 
lines on the appropriate graphic. The computed graph 
maximum and minimum are input from the graphic subroutine 
along with the axis tick intervals. Grid lines are output 
to the CRT from GRID, Lines 3900-4010 contain the 
subroutine. GRID requires 2 % of program space. 

This concludes the discussion of general 
graphics formating subroutines. The driving graphic 
subroutines will now be discussed. 

Subroutines GRIN, GRPOUR, and GRAM? are all 
similar in formatting structure. »ith some judicious choice 
of programming flags and routine definition, it should not 
be too difficult to combine these three subroutines into a 
single driving graphic subroutine. GROP is formatted 
differently, but still retains a few features common to all 
the graphic driving subroutines. The similarities between 
the subroutines will be discussed first. The particular 
differences applicable will be discussed secondly. 

Each graphic subroutine is composed of four 
basic blocks. They are: 

1. Input waveform data from mass btorage 

2. Define waveform and graphic boundaries 
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3. 



Draw the graph of the function using special defining 



functions when necessary 
4. Labling the graph with the correct title (s) 

GEiPH DRIVER inputs the function value PO which 
is used to identify the location of the data storage file in 
which the appropriate data rests. This single row array is 
usually returned in X. The waveform window and viewport 
parameters, bounding the output graphic display, are next 
determined. The maximum horizontal axis values are found by 
mutiplying the distance between a single data point (in 
UDU's) by the total number of data points in the array along 
the X-axis. For those waveforms listing the mean graph 
values as part of the descriptive parameters, the sum of all 
the array element values is obtained and then divided by the 
discrete number of array points: 

N 

Mean = W5 = E x(i)/N [4.11) 

i=l 

To ensure that the function is graphed along the zero axis 
of the formatted graphic, the means are subtracted from the 
array elements. In similar fashion, the root mean square is 
found: 
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( 4 . 12 ) 



RMS = W6 = [ Z x(i)^/512]^/^ 
i=l 



The maximum and minimum values for the input 
data are found by implementing the 4052 call '’MAX*' and "MIN'* 
routines. These values are returned in W4 and M3 
respectively. The values are then output to the axis 
defining subroutines for establishment of appropriately 
scaled and labled axis. 





Next , 


the GS cursor is positioned to 


the 


beginning 


of the array with a 


”M07E'’ command. For 


the 


number of 


discrete 


data in the 


array, the function 


is 


graphed on 


the GS C3 


T. Finally, 


the graphic is titled 


and 


identified 


by the 


appropriate 


parameters. Figure 4 


. 10 



provides a general flow diagram of the graphic driver 
routin es . 

The particular aspects of each graph will now be 
briefly discussed. 

f. Subroutine GRIU 

This subroutine is the driving graphic for the 
display of the inputted Direct, Incident and Augmented 
Waveforms, and derived Backscatter waveform. It resides in 
statements 5000-5410 and makes-up 7% of the program. 
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The general description of graphic routines 
provided above is most typified by GRIN. 

Two special notes need to be made. First, 

graphing begins with the ninth element of the array. This 
is to prevent the display of an erroneous spike input by the 
DPO. (These erroneous values have a minimal effect on the 
processing). Secondly, the function is graphed in discrete 
steps equal to l/H of the total array value along the axis: 

Graphed Point = (Q * 21), X (Q) CH. '1) 

where: Q = array number 

X (Q) = array element 

The title of the graph is selected by the 
routine located in lines 5290-5380. 

The horizontal dimenion for the graphic is time 
while the vertical axis is scaled in volts, 
g. Subroutine GRFOUR 

This subroutine is also fairly typical of the 
general graphic method. It displays the magnitude and phase 
components of the spectral data in polar fora for the Direct 
and SacXscattered waveforms and for the impulse response of 
the target. 
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Only the first 60 elements of the magnitude 
array are graphed. Values beyond the first 60 are of a 
random nature due to the processing performed to derive the 
transforms. Spectral content ouch beyond 10 GHZ is, 

therefore, negligible. All values of the phase array are 
plotted, although the note concerning the randomness of the 
data is also applicable here. 

GRFOUR is actually performed in two distinct 
operations. The first graphs the magnitude of the 

transform. The horizontal axis dimensions are converted to 
frequency from time by finding the frequency of each data 
point: 

Graphed Point = Q/(Z1*257), X (Q) 

The second operation graphs the transform phase. 
The horizontal axis dimensions are converted to radians: 

Graphed Point = (Q*27t)/257, X (Q) 

GSFOUR is programmed in lines 5500-6080. It is 10% of the 
total program. 

h. GROP 

GROP is the graphic that reproduces the shape 
representing the target. That is, the target shape based on 
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physical optics principals and uncorrectsd for field effects 



induced by residual currents. It represents a 
two-dimensional cross-sectional area cut of the target. 
Therefore, only target length and height are shown. 

This subroutine is constructed generally in the 
same format as the other graphic subroutines. One important 
difference is that a true one-to-one relationship must be 
esxablished for the coordinate axis. However, before this 
is done, the metric dimensions of the target length must be 
bounded. This is done in lines 6210-6240. Note that the 
array BO, source for the shapes storage in MATH, has already 
been scaled for the metric dimension in the y-direction or 
amplitude. The x-dimension, or length, is scaled for metric 
measures in line 6220. This is accomplished by multiplying 
the time between samples by the total number of samples in 
the array representing the target. This value is then 
compensated by the averaging antenna parameter established 
in MATH, S5: 



- (2)(L)(1.2xlO^°HSS') 

Q (4.1 

where: Z = sampling time between single samples 

L = measured length in number of samples 
S5 = antenna length compensator 
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Q = low-pass filter compensator 
Q is always one for zero optimization, or no low-pass 
filtering. It has effect only when optimization is 

performed, and is designed to compensate for the spreading 
of the target dimensions in the x-direction during 
filtering. Note also, that all defined targets are measured 
in inches. All outputs are scaled inches. It should be a 
relatively simple procedure to generalize the program to 
handle either metric or english units. 

The next major difference between GSOP and other 
graphic routines is that all graphing is done in GROP itself 
without calling AXIS or any other external formatting 
subroutine. As mentioned above, the important concept is to 
properly relate the co-ordinant axis on a true one-to-one 
plot. This is done is statements 6290-6360 and is called 
•’axis scaling”. Axis scaling is done in two steps. First 
the general x-axis scaling factor AO is found. This is done 
in 6290 by dividing the maximum horizontal bounding UDU by 
the viewport established by the display; 

AO = (W2 - W1) / (V2-V1-5) 

This yields DDU's/GDO in the x-direction. The y-axis 
scaling factor A1 is then fixed to this value. The 
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dimensions in the y-direction are then checked by the 
program in lines 6320 to ensure that ’’clipping” or windowing 
out of the data does not occur. This is necessary as the 
vertical viewing axis is only 62H that of the x-axis. If the 
y-dimension were more than S2% of the x-dimension, data 
would be lost during graphing. To prevent this, the x-axis 
values are fixed so that the maximum x-dimension will always 
be plotted. This is an effective reduction in the size of 
the total displayed shape. To retain the true one-to-one 
axis relationship, the y-direction dimensions must be 
reduced by a corresponding percentage. This reduction 
factor is computed in line 6330 and applied to the array in 
6440. Note that in 6330 twice the maximum y-dimension is 
used. This is done because the graphic output produces the 
mirror image of a axisymmetric object while the array itself 
is single-sided. 

Once scaling is established in statement 6380, 
the array is plotted. As data is graphed it must be 
dimensioned in english units, compensated for optimization 
effects, and compensated for any viewport clipping. This is 
done on a ?oint-by-point basis in line 6410. Note that the 
array has already been totally dimensioned and compensated. 
It is in fact the plotting point which now must be taken 
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into consideration. This compensated plotting point is 
returned in 01 and is used in drawing the graph. Finally, 
the mirror image of each point is graphed to produce a 
"true" image. A flow diagram for GROP is given in Figure 
4. 1 1. 

i. Subroutine GRAHP 

GRAMP graphs the time domain ramp response in 
lines 6600-6910. It comprises 5J5 of the program. 

GRAMP is identical to GRIN in practically every 
detail. Note that the array X is multiplied by 1000 in line 
6670. This converts array amplitudes to millivolts. 
GRIN also employs this conversion. 
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In this chapter, the results of several selected 
measurecients vill be presented. The observations made were 
chosen to demonstrate the flexibility of the Operating 
System as well as verify experimental results with the 
theoretical. 

k. aEASOHEHEHT 1 — 7ABT BOHESISHTED TARGET DISTANCE 
♦ Descr iption of Setup 

The purpose of this measurement is to observe what 
effect locating the target a distance different from the 
transmitting antenna then that specified for the receiving 
antenna, but maintaning the boresighted relationship of the 
antennas has on measurements taken by the OS. The 

receiving antenna is placed 2.54 l! from the transmitting 
antenna. The target is defined from the target library on 
Tape 2. 

2 . Result s 

The results are given in Figures 5.1 to 5.3 with the 
backscattered waveform, ramp response and physical optics 
shapes being presented. 
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The ramp response is seen to be typical for a 
cylinder although somewhat reduced in size. This is to be 
expected since the propagation attenuation losses are 
greater for the distance at which the target is located in 
this experiment than it is for the normal control case. 
With the distance taken into consideration by the Operating 
System software# it can be seen that the resultant physical 
optics shape is a close approximation. The target width 
(length) is within 10% of the true length and the diameter 
is 15.4% of the true diameter. The uncorrected width is 49% 
greater than the true length. The backscattered waveform 
shows the effect of some spreading in comparison to the 
control case and also displays a reduced S/N. 

3 . C one lu si ons 

The distance between the receiving antenna and 
target can be arbitrary and the Operating System will be 
able to provide a reasonably accurate interpretation of the 
target. As the distance increases, this accuracy can be 
expected to decrease directly to the limit of the Systems 
ability to seperate the target and noise. 
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target: run date dist ant tgt remarks 

2 3-7-81 2.54N I 1 LOCATION DIFFERENCE 

NAXIMUN PEAK UALUE nU 

HININUn PEAK UALUE -1.48 nU 

RMS UALUE nU 

KEAN UALUE 9.09 nU 

NUMBER OF HAUEFORMS AUERAGED ■ 50 OPTIMIZATION UALUE - 8 

BACK SCATTER MAUE FORM 




Figure 5.1. MZasUREHENT 1 — Backscattered Wavefora. 
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TARGET: RUH DATE DIST AHT 

2 8-7-01 2.34H 1 

MAXinUH PEAK UALUE 

(1IHIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

NUMBER OF NAUEFORMS AUERAGED - SB 



TGT REMARKS 

1 LOCATION DIFFERENCE 
+6.86 fiU-M 
-196.33 nU-li 
♦112.76 nU-M 
-90.04 nU-M 

OPTIMIZATION UALUE - 9 




8.30 1.00 1.30 2.00 2.30 3.00 3.38 4.00 4.58 3.38 nS 



Figure 5.2. HEASOREHENT 1— Ramp Response 
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TARGET: RUN DATE DIST 

2 8-7-81 2.34M 1 

DIAMETER < IHCHES-CORRECTED) 
HIDTH <1HCHES-C0RRECTED) 

WIDTH <1HCHES-UHC0RRECTED)«* 
HUMBER OF WAUEFORMS AUERACED • SB 



ANT TGT REMARKS 
1 1 LOCATION DIFFERENCE 

8.9825 < 8.0000) 

13.2047 <12.0000) 

17.3711 <12.0000) 

OPTIMIZATION UALUE ■ 0 



RAMP RESPONSE ~ PHYSICAL OPTICS 




Figure 5.3. HEiSOREMENr 1 -—Physical Optics Shape 
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B. MEASOREHENT 2 — TARGET AND RECEIVING ANTENNA NOT 

BORESIGHTED 

1 . Description of Setup 

An attempt is made to observe the effect of placing 
a target at an arbitrary location on the IDRL ground plane. 
The transmitted field is not axially incident on the target. 

For this experiment, the target (cylinder 1) was 
placed 2.11 meters from the transmitting antenna and at a 
90® angle (approximately) to the recaiving/transmirting 
antenna axis. The receiving antenna is located in the 
control area of the TDRL. The target is at about a 60® 
angle to the receiving antenna and located 2.44 meters from 
it. The experimental setup is diagramed in Figure 5.4. 

2. Re su It s 

The results are given in Figures 5.5 to 5.7 
representing the ramp response, physical optics shape and 
backscattered waveform. Note that the distance input to the 
system is incorrect. The correct value as read from the DPO 
should be 1.57 meters. The diameter should be multiplied by 
the factor 1.53. Doing so yields a diameter of 6.3341 
inches as the true corrected value. The resulting 

compensated values are thus very close to the targets actual 
dimensions. These figures may be suspect as subsequent runs 
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Figure 5.4. Measureaeat 2--Physical Setup 

never yielded the saae degree of accuracy, although they 
approached it. The important item is that the physical 
optics shape is a close approximation of that expected to be 
representative of a cylinder. 

3 • Conclusion s 

It would appear that the targets location can be 
generalized if the shape only and not the dimensions is the 
significant item of interest. However, time precluded 
further verification of this observation. It cannot be said 
with absolute certainty that this observation is completely 
general. other configurations with the target located at 
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different distances from the boresighted transmitting and 
receiving antennas should be observed. Note that for proper 
compensation for propagation attenuation, the distance the 
target is from the receiving antenna is determined by 
reading the DPO Distance Scale. This should be the OS 
distance input value vice the straight line measured 
distance from transmitting antenna to target. 
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target: run date dist ant tgt 

3 8-7-91 1.27H 1 I 



REMARKS 

NOT BORESICHTEO 



MAXIMUM PEAK UALUE +47.50 mU-M 

MINIMUM PEAK VALUE -329.05 nV-M 

RMS UALUE +181.31 «U-M 

MEAN UALUE -133.00 «V-M 

NUMBER OF NAUEFORMS AVERAGED • 4 OPTIMIZATION UALUE 




0.50 1.80 1.58 2.00 2.50 3.00 3.50 4.00 4.50 5.00 nS 



Figure 5.5. 



Heasureneat 2— Samp Response 
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TfiO^FT* f^LlM DATE DIST ANT TCT REMARKS ■' **t*^i* 

TARGET. RUN DATE^^ ^ ^ BORESIGHTED 

DIAMETER < INCHES-CORRECTED) ,t*A?l9 fl2*0000> 

iiiVn I is'csiiiss^si i’ . irs^LUE . , 

HUMBER OF HAUEFORMS AUERAGED^»_4 

response — PHYSICAL OPTICS 




Figure 5.6 



Measurement 2 — Physical Optics Shape 



target; run date D13T AHT TCT 

3 8-7-01 1.27M 1 I 



REMARKS 

HOT BORESICHTEO 



MAXIMUM PEAK UALUE 

MINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

HUMBER OF WAUEFORMS AUERAGED - 4 



•►2.48 nU 
-2.49 nU 
•►0.98 nU 
8.00 nU 

OPTIMIZATION UALUE 



BACK SCATTER HAUE FORM 



nU 

2.80 
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Figure 5.7. aeasureaent 2 — BacJcscattered flavefora 
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C. aZaSOEBMENT 3 — AHTEHNA EFFECTS OH TARGET RESPOHSES 

1 . De scri ptio n of Se tup 

T5te purpose of this measurement is to verify that 
the effects of the antenna are in fact deconvolved out 
during processing. aost previous runs have been performed 
using antenna 1 (TEH Horn) . This measurement is performed 
using three different length monopole antennas. The 
antennas used were antenna 2 (see TABLE II) (85 mm 

aonopole) , antenna 3 (205 mm monopole) and antenna 12 (285 
mm monopole) . The target used was target 1 (see TABLE I) . 
Both target and antenna were placed in the control 
locations, 1.27 meters from the transmitting antenna. 

21 -acquisitions were averaged for each run. The results 
obtained should be compared with those obtained previously 
for rhe TEH horn. 

2 . Resu It s 

The results are presented in Figures 5.7 to 5.24 
Each antenna has the results of its run shown in six graphs. 
The first graph depicts the direct waveform as seen by the 
individual antenna, the second is the ramp response 
resulting, the third figure is the physical optics 

representation of the target shape, and the final three 
figures are the optimized frequency domain impulse, time 
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,:s: 












i^'-'V 







domain impulse and ramp response for the target. Figure 5.7 
to 5.12 are for monopole antenna 2. Figures 5.13 to 5.18 
are for monopole antenna 8 and figures 5.19 to 5.24 are for 
aonopole antenna 12. 

The first obvious result is that in comparison of 
the ramp responses, the smaller rhe antenna is in relation 
to the target, the less typical is the ramp and, 
correspondingly, the physical optics shape. However, for 
antennas considerably taller than the height of the target, 
a significant amount of noise pollutes the ramp response. 
In comparing the target impulse responses, it can be seen 
that the exact opposite appears to be true. The most 
representative impulse is for the the shortest monopole. 
The frequency domain impulse graphs show that the longer 
monopoles accept a higer frequency content than the shortest 
monopole. 

3 . C onclusions 

First, the closer the antenna is to natural 
resonance with the incident pulse, the more faithful appears 
the resulting impulse response to theory. But it does not 
naturally follow that the corresponding ramp response will 
be as true a r epsentat ion. The opposite appears to be true. 
It may be that the short antenna is shadowed by a portion of 
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th-e target. More likely, the short antenna is just not 
large enough to accept the full radiated backscatter of the 
target. Conversely, the 235-mm monopole is sensitive to 
environmental noise. As can be seen in figures 5.20 and 
5.14, the ramp responses display the relatively greater 
sensitivity of the 285-mm monopole to environmental noise 
than the 205-mm monopole. Although the observation is 
limited to just these runs, it may be possibly concluded 
that there is a "best" monopole antenna length for a 
particular target size. 

In conclusion, it can be said that the effects of 
the antenna are in general deconvolved out. However, other 
effects indirectly related to the target/antenna 
configuration may introduce factors which materially effect 
the reproduction of the target shape and its overall 
responses. 



175 



TARGET; RUN DATE DI3T AHT TGT' REMARKS 

13 8-11-81 1.27M 2 1 DIPOLE ANTENNA #2<85MH) 

MAXIMUM PEAK UALUE +76.01 nU 

MINIMUM PEAK UALUE -91.33 nU 

RMS UALUE ♦SI. 33 hU 

MEAN UALUE 0.00 mU 

NUMBER OF WAUEFOPMS AUERAGED =• 21 OPTIMIZATION UALUE « 0 




0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 nS 



Figure 5.8. Heasareaent 3 — Direct Waveform (85 am 



176 



TARGET: RUN DATE DIST 
13 3-11-81 1.27M 

tIAXIMUH PEAK UALUE 

MIHinUN PEAK UAUUE 

RMS UALUE 

MEAN UALUE 

NUMBER OF HAUEFORNS AUERAGED 



ANT TGT REMARKS 

2 1 DIPOLE ANTENNA i2<85NN> 

.. <-782. 43 nU-M 

.. -503.93 mU-M 

♦374.80 mU-N 
.. ♦163.01 nU-M 

21 OPTIMIZATION UALUE « 0 




8.50 1.88 1.50 2.88 2.50 3.88 3.50 4.88 4.58 5.80 nS 



Figure 5.9. Heasureaent 3 — Ramp Response (85 mm) 



target: ruh date oist 

13 8-11-91 1.27M 2 

DIAMETER < IHCHES-CORRECTEO> 
WIDTH (INCHES-CORRECTED^ 

WIDTH <1HCHES-UNC0RRECTED)»» 
HUMBER OF WAUEF0RM3 AUERAGED - 21 



ANT TGT REMARKS 

2 1 DIPOLE ANTENNA #2<85MM> 

5.5879 < 5.8909) 

12.3283 <12.0000> 

10.1367 C12.0000) 

OPTIMIZATION UALUE ■ 9 



RAMP RESPONSE -- PHYSICAL OPTICS 




Figure 5,10 



aeasure*ent 3--Ph.ysical Optics shape (85 mn) 



TARGET: RUN DATE DIST ANT TGT REMARKS 

13 8-11-91 1.27M 2 1 DIPOLE ANTENNA #2<85MI1) 

MAXIMUH PEAK UALUE 8.0008 nU 

MINIMUM PEAK UALUE -0.0007 nU 

RMS UALUE 0.0002 nU 

MEAN UALUE 0.0000 nU 

HUMBER OF WAUEFORMS AUERAGED - 21 OPTIMIZATION UALUE « 0.3 



nU 

0.0005 
0.0000 
0.0005 

0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 nS 



IMPULSE RESPONSE 




Figure 5.11. 



Measuremenl: 3 — Time-Domain Impulse Response 
(85 mm) 
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TARGET: RUN DATE DIST ANT TGT REMARKS 

13 8-11-81 1.27M 2 1 OIPOLE ANTENNA #2<85MM) 

NUMBER OF WAUEFORMS AUERAGEO » 21 OPTIMIZATION UALUE ■ 0.5 



IMPULSE FFT— MAGNITUDE 





0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 E.0 RAO 



Figure 5.12. 



Measurement 3 — Preguency-Domain Impulse Res- 
ponse (8 5 mm) 
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TARGET: RUN DATE DIST 
13 3-11-81 1.27M 

MAXIMUM PEAK UALUE 

MINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

NUMBER OF HAUEFORMS AUERAGED 



ANT TCT REMARKS 

Z 1 OIPOLE ANTENNA «2<85MM> 

+613.71 nU-M 

. ... -345.20 nU-M 

.... +284.65 nU-H 

. . . . +141.70 nU-M 

» 21 OPTIMIZATION UALUE * 0.5 



RAMP RESPONSE — TIME 




Figure 5.13. Seasureaent 3 — Ramp Optimized (85 ma) 



TftRGET: RUH DATE DIST ANT TGT REMARKS 

14 8-11-81 1.27M 8 1 DIPOLE ANTENNA 18 <205MH) 

MA.KIMUM PEAK VALUE +104.92 mU 

MINIMUM PEAK VALUE -116.61 (iU 

RMS VALUE +59.50 nV 

MEAN VALUE 0.00 tiU 

NUMBER OP WAVEFORMS AVERAGED » 21 OPTIMIZATION VALUE * 0 




8.58 1.00 1.58 2.80 2.50 8.08 3.50 4.00 4.50 5.88 nS 



Figure 5, 14. Measureaent 3 — Direct Waveform (205 am) 
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ThiRGET: RUN DATE 

14 8-11-91 



DIST ANT TCT 
1.27M 8 1 



DIPOLE ANTENNA «8 <205MN) 



MAXIMUM PEAK UALUE +5.33 nU-M 

MINIMUM PEAK UALUE -609.89 nU-M 

RMS UALUE +357.09 nU-M 

MEAN UALUE -299.54 mU-M 

NUMBER OF WAUEFORMS AUERAGED » 21 OPTIMIZATION UALUE 




8.58 1.80 1.58 2.80 2.58 3.80 3.50 4.30 4.50 5.00 nS 



Figure 5.15. 



Measurement 3— Ramp Response (205 mm) 
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target: run date oist 

14 8-11-81 1.27M 

DIAMETER ( INCHES-CORRECTED) 
WIDTH (INCHES-CORRECTED) 
WIDTH <INCHE3-UHC0RRECTED)*» 
HUMBER OF WAVEFORMS AVERAGED » 



AHT TGT REMARKS 

8 I DIPOLE AHTEHNA #8 <2051111? 
3.3341 < 6.0000) 

12.0734 <12.0000) 

19.3945 <12.0000) 

21 OPTIMIZATION VALUE - 0 



RAMP RESPONSE — PHYSICAL OPTICS 




Figure 5.16 



Seasureaent 3 — physical Optics Shape (205 am) 



TARGET; RUH DATE DIST AMT TGT REMARKS 

14 8-11-91 1.27M 8 1 DIPOLE AMTENHA «S (203MI1) 

MAXIMUM PEAK UALUE 8.0005 nU 

MINIMUM PEAK UALUE -0.0006 nU 

RHS UALUE 0.0002 nU 

MEAN UALUE 0.0000 hU 

NUMBER OF WAUEFORMS AUERAGED » 21 OPTIMIZATION UALUE - 0.5 



nU 

0.0005 
0.0000 
-0.0005 

0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 nS 




Figure 5.17, 



Heasureient 3 — Time-Domain Impulse Response 
(205 mm) 
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TARGET: RUN DATE DIST ANT TGT REMARKS 

14 8-11-81 1.27M 8 1 DIPOLE ANTENNA 18 <205M«) 

HUMBER OF WAUEFORMS AUERAGED « 21 OPTIMISATION UALUE - 0.5 



IMPULSE FFT— MAGNITUDE 





0.8 8.5 1.0 1.5 2.8 2.5 3.8 3.5 4.8 4.5 5.8 5.5 6.0 RAD 



Figure 5.18. fleasureaent 3 — Frequency-Domain Impulse Res- 
ponse (2 05 mm) 
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target: ruh date dist ant 

14 8-11-31 1.27M 8 

MAXIMUM PEAK UALUE 

MINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

NUMBER OF WAUEFORMS AUERAGED * 21 



tct remarks 

1 DIPOLE ANTENNA 18 <283MH) 
♦0. 19 mU-M 
-537.10 fiU-M 
♦318.09 mU-M 
-255.31 mU-M 

OPTIMIZATION UALUE » 0.5 




0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.09 nS 



Figure 5. 19. Measurement 3 — Ramp Optimized (205 mm) 
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TARGET: RUN DATE DIST 
12 8-11-81 1.27M 

riAXiriUM PEAK UALUE 

tllHIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

HUMBER OF WAUEFORMS AUERAGED » 21 



REMARKS 

DIPOLE ANTENNA •12(285Mri> 
♦101.31 nU 
-122.78 fiU 
♦67.32 nU 
0.00 nU 

OPTIMIZATION UALUE * 0 



ANT TGT 
12 1 




0.50 1.88 1.50 2.00 2.50 3.00 3.50 4.00 4.58 5.80 nS 



Figure 5.20. Measurement 3 — Direct Waveform (285 mm) 
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target: run date oist ant 

12 8-11-31 1.27M 12 

MAXIMUM PEAK UALUE 

MINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

HUMBER OF WAUEFORMS AUERAGED » 21 



TCT REMARKS 

1 DIPOLE ANTENNA »12<28SMM> 

+ 139.35 tiU-M 

-563. 13 nU-M 

+329.92 mU-M 

-252.72 mU-H 

OPTIMIZATION UALUE * 9 




9. 59 1.98 1.59 2.09 2.59 3.99 3.59 4.98 4.59 5.99 nS 



Figure 5 . 21 . Heasurenent 3 — Samp Response (285 mm) 
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TARGET; RUN DATE DIST ANT TGT REMARKS 

12 8-11-91 1.27M 12 1 DIPOLE ANTENNA #12<293MH) 

DIAMETER < INCHES-CORRECTED) €.8000 ( 6.9000> 

WIDTH UNCHES-CORRECTED' 12.8000 <12.0000> 

WIDTH <INCHES-UNCORRECTED)t» 15.2303 <12.0000> 

NUMBER OF WAUEFORMS AUERACED » 21 OPTIMIZATION UALUE ■ 0 



RAMP RESPONSE — PHYSICAL OPTICS 




Figure 5.22. Measuceaent 3 — Physical Optics Shape (285 mm) 



TARGET: RUN DATE C>IST ANT TGT REMARKS 

12 8-11-81 1.27M 12 I DIPOLE ANTEHHA *12<285HM) 

MAXIMUH PEAK UALUE 0.0006 mU 

MINIMUH PEAK UALUE -0.0003 nO 

RMS UALUE 0.0002 nU 

MEAN UALUE 0.0000 mU 

NUMBER OF WAUEFORMS AUERAGED » 21 OPTIMIZATION UALUE « 0.5 



nU 



0.0005 



0.0000 



0.0005 




0.50 1.00 1.50 2.00 2.50 3.00 3.30 4.00 4.58 5.00 nS 



Figure 5.23 



Measureaent 3— Tiae-Donain lapulse Response 
(285 am) 



TARGET: RUN PATE PIST AHT TGT REMARKS 

12 3-11-81 1.27M 12 1 DIPOLE ANTENNA I12(285HM> 

NUMBER OF WAUEFORMS AUERAGED » 21 OPTIMIZATION VALUE • 0.3 



IMPULSE FFT— MAGNITUDE 





0.0 0.3 1.0 1.3 2.0 2.3 3.0 3.3 4.0 4.3 3.0 3.3 £.0 RAD 



Figure 5.24. Measurement 3 — Frequency-Domain Impulse Res 
ponse (2 85 mm) 
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TARGET; RUN DATE DIST ANT 
12 8-11-91 1.27M 12 

MAXIHUn PEAK UALUE 

MIHIMUH PEAK UALUE 

RMS UALUE 

MEAN UALUE 

NUMBER OF WAUEFORMS AUERAGED ■ 21 



TGT REMARKS 

1 DIPOLE ANTENNA il2<285MM> 
+28.80 hU-M 
-59S.67 nU-M 
+368.85 mU-M 
-314. 14 tiU-t1 

OPTIMIZATION UALUE - 0.5 




0.S0 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 nS 



Figure 5.25. aeasurement 3--Raaip Optimized (285 mm) 



D. MEASORESENT 4 — ARBITRARY COMPOSITE TARGET 



1 . Descript ion of Setup 

This measurement was made to observe the response of 
the OS to a completely arbitrary, but axisymmetric target 
not encountered by the system previously. The composite 
target was a combination of target 11, a cone, and target 3, 
a cylinder. A drawing of the composite is given in Figure 
5. 25. 




Figure 5.26. Measureaent 4 — Composite Target 
2 . Re suit s 

The results are given in Figures 5.26 and 5.27, the 
physical optics shape and ramp responses respectively. The 



corrected diameter 


is 


11 % 


larger than actual, and 


the 


corrected length 


is 


35% 


shorter than actual. 


The 
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uncorrected length is 13 % greater than actual. The shape is 
a close approximation of the targets actual shape, although 
somewhat shortened. The cone portion appears to be about 
one-half of the the total figure, as it was in the actual 
target . 

3 . C onclusion s 

Based on this observation and the many others taken 
during the implementation of the OS, it can be said that the 
OS can be relyed upon to produce recognizable representative 
shapes for known or unknown axisymmetric metallic targets. 



195 



TARGET: RUN DATE DIST ANT TGT REMARKS 

27 9-11-31 1.27N 1 ??? COMPOSITE ARBITRARY TARGET 

DIAMETER < INCHES-CORRECTEO) 6.5639 ( 0.0000) 

UIDTH (IHCHES-CORRECTED) 8.9094 ( 0.0000) 

WIDTH CINCHE3-UNC0RRECTED)t< 13.5352 ( 0.0000) 

UMBER OF HAUEFORMS AUERACED » 21 OPTIMIZATION UALUE - 0 



RAMP RESPONSE — PHYSICAL OPTICS 




Figure 5.27 



Heasuremeut 4~—Physical Optics Shape 



target: run date dist 

27 3-11-31 1.27H 

MAXIMUM PEAK UALUE 

MINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

NUMBER OF WAUEFORHS AUERAGED 



ANT TCT REMARKS 

1 7?? COMPOSITE ARBITRARV TARGET 



■*•38.86 nU-M 
-335.77 nU-M 
+171.95 mU-M 




-115.28 nU-M 



21 OPTIMIZATION UALUE » 8 




0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5. 00 ti$ 



Figure 5.28. Measurement 4 — Ramp Response 
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E. MEASDREHENT 5 — DIELECTBIC CTLIHDERS 



1 . Descript ion of Setup 

Two runs were aade using different dielectric 
cylinders in place of metallic objects as targets. The 
dielectric cylinders used are similar in configuration to 
cylinder targets 2 and 3 machined from aluminum stock, but 
are made of plastic. For this measurement, the dielectric 
cylinders will refered to as target 2 and 3, respectively. 

The purpose of the measurement was to observe the 
response of the dielectric cylinders to axially incident 
electromagnetic fields and to also see how the OS responds 
to a target which was non-metallic. The targets were placed 
in the standard control location, 1.27 meters from the 
transmitting antenna. The runs assumed the targets were of 
arbitrary, unknown shape for DS observation. 

2 . R esult s 

Target 2 was run first. The results are given in 
Figures 5.28 to 5.31. The physical optics shape in Figure 
5.23 is a very poor representation of the true targets 
shape. The reason for this can be observed in the ramp 
response of Figure 5.29. The S/N is extremely low. A noise 
spike at about 1.6 nsecs has caused the entire resultant 
shape to be foreshortened. The amplitude discrepancy can be 
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attributed to the fact that the amplitude compensator was 
developed for the relatively large valued returns of the 
metallic cylinders. Nevertheless, when the ramp response is 
optimized, the general shape of the target can be seen as in 
Figure 5.30 and 5.31. 

The results using Target 3 are shown in Figures 5.32 
to 5.35. Because this target is physically larger than 
target 2, it has returned a substantially larger signal, 
large enough so that system noise is less of a significant 
factor (although still seen on the waveform) . The amplitude 
is again much below that to be expected for a metallic 
target of similiar size. 

3 • Co nclus ions 

Although the information provided by this 
measurement is limited, it appears that dielectric targets 
might provide returns which can be identified by time domain 
techniques. This could be an important area for more 
intensive research as the use of high density polymers and 
resins in real targets of interest is of growing importance 
to surveillance systems attempting to locate and identify 
such targets. 
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target: run date oist ant tct remarks 

24 9-11-31 1,27(1 1 ??? DIELECTRIC CYLINDER <3X12 IH) 

DIAMETER ( INCHES-CORRECTED) 0.7662 < 8.0000) 

WIDTH (INCHES-CORRECTED) 2.1596 < 0.0000) 

WIDTH <INCHES-UNCORRECTEDit» 3.2813 < 0.0000) 

NUMBER OF WAUEFORMS AUERAGED ■ 21 OPTIMIZATION UALUE - 0 

RAMP RESPONSE — PHYSICAL OPTICS 




Figure 5.29 



Measureaent 5--Physical Optics (Target 2) 



TARGET: RUN DATE DIST ANT 
24 8-11-31 1.27H 1 

MAXIMUM PEAK UALUE 

MINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

HUMBER OF WAUEFORMS AUEPAGED - 21 



wU-M 

10 

5 

0 



-10 

-15 

-20 

-25 

-30 

-35 



TCT REMARKS 

??? DIELECTRIC CYLINDER <3X12 IN) 
>11.15 mU-M 
-38.59 fiU-M 
+17.07 mU-M 
-13.94 mU-M 

OPTIMIZATION UALUE « 0 



RAMP RESPONSE — TIME 




0.50 1.00 US0 2.00 2.50 3.00 3.50 4.00 4.50 5.00 nS 



Figure 5.30. Heasarement 5--Ramp Respoase (Target 2) 
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target: run date diet ant tgt remarks 

24 9-11-81 1.27M 1 ??? DIELECTRIC CYLINDER <3X12 IN) 

DIAMETER ( INCHES-CORRECTED) 0.7662 < 0.0000) 

WIDTH < INCHES-CORRECTED) 2.1596 < 0.0000) 

WIDTH (INCHES-UNCORRECTED)tt 9.0820 < 0.0000) 

NUMBER OF WAUEFORMS AUERAGED • 21 OPTIMIZATION UALUE - 1 



RAMP RESPONSE — PHYSICAL OPTICS 




Figure 5.31 



Measureaent 5“~3ptiinized Physical Optics 
get 2) 



(Tar- 



202 



TARGET: RUN DATE DIST ANT 
24 3-11-81 1.27N 1 

MAXinUM PEAK MALUE 

MINIMUn PEAK UALUE 

RHS UALUE 

MEAN UALUE 

NUMBER OF WAUEFORHS AUERAGED = 21 



TGT REMARKS 

??? DIELECTRIC CYLINDER <3X12 IH> 
+8.29 nU-M 
-27.71 liU-N 
+13.28 nU-M 
-10.39 mU-H 

OPTIMIZATION UALUE « 1 ' 



RAMP RESPONSE — TIME 







Figure 5.32, Measurenent 5--Optimized Raap (Target 2) 
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TARGET: RUN DATE DIST ANT TGT REMARKS 

23 S-H-31 1.27M 1 ??? DIELECTRIC CYLINDER <6X6 IH> 

DIAMETER < INCHES-CORRECTED) 1.4134 < 0.0000) 

WIDTH «; INCHES-CORRECTED) 5.7075 < O.0000) 

WIDTH <INCH£3-UNC0RRECTED>t* 3.6719 < 0.0000) 

NUMBER OP WAUEFORMS AUERAGED - 21 OPTIMIZATION UALUE - 0 



RAMP RESPONSE — PHYSICAL OPTICS 




Figure 5.33. 



ileasur ement 5--Physical Optics (Target 3) 
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target; run date oist ant 

25 8-11-81 1.27M 1 

MAXinUN PEAK UALUE 

MINIMUM PEAK UALUE 

RMS UALUE 

MEAN UALUE 

NUMBER OF WAUEFORMS AUERAGED « 21 



TGT REMARKS 

DIELECTRIC CYLINDER <6X6 IH> 
♦94.83 mU-M 
-71.29 mU-M 
♦35.00 mU-M 
-3.50 mU-M 

OPTIMIZATION UALUE • 9 



RAMP RESPONSE — TIME 




Figure 5.34. Measarement 5“-Ramp Response (Target 3) 
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target: run date dist ant tgt remarks 

25 8-11-31 1.2TM 1 ??7 DIELECTRIC CYLINDER <6X6 IN) 

MAXIMUM PEAK UALUE 8.0983 nU 

MINIMUM PEAK UALUE -8.0884 hU 

RMS UALUE 8.0881 nU 

MEAN UALUE 0.0000 nU 

NUMBER OF WAUEFORHS AUERAGED » 21 OPTIMIZATION UALUE ■ 9.5 



nU 



8.8882 



0.8888 



-8.8882 




8.58 1.88 1.58 2.00 2.50 3.90 3.50 4.08 4.58 5.98 nS 



Figure 5.35. Measureaent 5 — Dptimized Physical Optics (Tar- 
get 3) 
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TftRCET: RUN DATE OIST AMT TCT REMARKS 

25 8-11-81 1.2711 1 ??? DIELECTRIC CYLIHDER <6X6 IH> 

NUMBER OF WAUEFORHS AUERACED ■ 21 OPTIMIZATION UALUE - 0.5 



IMPULSE FFT— MAGNITUDE 





0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4,5 5.0 5.5 6.0 RAD 






Figure 5. 36 . 



Measarement 5--Optimized 8aap (Target 3) 
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SUMMARY 



A software Operating System has been developed for the 
Time Domain Radar Laboratory. The System has been designed 
to be relatively easy for inexperienced operators to use 
after an initial introduction zo the OS code and rhe 
constituent hardware controlled. The OS extends the 
capability of the TDRL to the point where accurately 
processed transient time-domain data is presented in a form 
suitable for operator evaluation. The OS should serve as a 
primay research tool for the srudy of electromagnetic 
transients at the Naval Postgraduate School. 

The OS is highly opera ror- processor interactive. 
Although the code that makes-up the OS is relatively 
complicated, irs actual implementation and use is not 
difficult. There are numerous directives issued by the 
System and displayed on the Graphic Systems Crt to the 
operator requesting interactively the input of information 
through the GS keyboard at vital points in the programs. 

The OS is composed of three subprograms. They are: 

1. INPUT for the inputting of signal information from 
the TDRL ground plane 



208 



2. MATH for psrforaiing the significaat mathematical 
computations for the OS 

3. GRAPH for the presentation of the processed 
information in a form suitable for user evaluation 
and analysis 

INPUT receives the desired signal in discrete digitized 
form from the systems Digital Processing Oscilloscope, the 
principal analog to digital converter used in the TDRL. The 
source signal is a rapid rise time gaussian pulse simulating 
an impulse. The acguired information is averaged a number 
of times as interactively requested by the operator. This 
is done to reduce the white gaussian noise content of the 
acquired analog signal. 

HATH performs the Fourier conversion and signal 
deconvolution to determine the impulse response of the 
target. It also implements a noise reduction algorithm to 
optimize the processed data. The amount of optimization is 
under the control of the operator and is dependant on the 
degree to which loss of small surface features can be 
tolerated and yet obtain a good representation of the 
targets true shape and dimensions. The step and ramp 
responses are also obtained and stored for later evaluation. 
HATH determines the Physical Optics Shape of the target 
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bas- 2 d on the ramp ras 
induced currants on t 
computes and applies 
width based on hi 
pre-defined targets, 
to unknown targets ras 
in determination of th 
GRAPH outputs the 
4052 Graphic System, 
time direct, incid 
waveforms, the time-d 
shape of the target a 
for the direct and ba 
frequency domain impu 
the optimization vari 
control to INPUT to ac 
The total result 
transient analysis to 
internal organization 
modification so as t 
hardware and softwar 
aspects of time-doma 
demonstrated in the co 



ponse uncorrected for the effects of 
he body of the targer. It further 
compensation factors for length and 
storically observed results from 
The compensators are applied directly 
ulting in an accuracy of 10% or better 
e true target dimensions. 

processed data to the GS CRT of the 
Information graphed includes the real 
ent, augmented and backsca ttered 
cmain ramp response, physical optics 
nd the frequency domain spectral data 
ckscattered waveforms as well as the 
Ise response. A means for inputting 
able and for the return of program 
quire a new signal is also provided. 

is a highly versatile and powerful 
ol. The OS has flsxiblity in its 
to allow for future expansion and/or 
o allow the control of constituent 
e in applications involving other 
in target identification. This is 
mplete integration of a Frony's Method 
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Hybrid Solution for the location of a 



targets poles 



and 



residues with the requirements for pure transient 

time-domaia analysis procedures. 

The code itself is complete as it stands. However/ there 
are areas in which it could be made mere efficient. 
Particularly/ a number of graphic subroutines could be 
combined into one routine so as to reduce memory 
requirements. klso, a method for the display and selection 
of menu items using placemens of the screen cursor to select 
the item might be developed. This would allow for the 
display of more informarion graphing processes than is new 
allowed. 

One of the specialized algorithms developed for the OS 
is a method for expanding certain displayed graphs while 
retaining a true plot of relevant information. The result 
is that the GS has essentially infinite resolution and would 
allow for a close look at particular areas cf the display 
without disturbing the true graph relationships. Although 
presently applied to only one subrourine in GRAPH/ it could 
easily be adapted to any other graphic subroutine. 

A second specialized routine is one whereby a true 
one--o-one relationship between coordinate axis is 
maintained. The unique feature of this algorithm is its 
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generality in being applicable no any size iana field on any 
size viewport without ever losing inforaation due to 
"viewport clipping". 

The results of field observanions using "^he OS have 
proven remarkably faithful to theoretical results. However, 
the present use of the OS has been principally directed 
towards verifying that the sysnetn was up and working. As 
such, data was collected for axisymmetric targets almost 
exclusively as there was a large amount of literary support 
in this research area. There srill is a need to explore 
system operation with this configuration of targets. In 
particular, the OS needs to be extended so as to be aola to 
solve the Magnetic Field Integral Equation for EM transients 
on the surface of a target. This calculation would then be 
compared with field results to obtain a corrected and 
accurate representation of the target both in shape and 
di mension. 

Hardware has proven very well suited to the needs of the 
TDRL. However, a monostatic antenna more directional than 
the present omnidirectional radiator presently used, may 
make the system more generally applicable. Improvements in 
reception by better match of the systems antenna components 
should allow for the determination of finer feature on 
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smaller targets. A mors rapid rise time generator having 
greater power outpats than prssently supplied by the UH? 
impulse generator would be useful. Finally, a fully 
computerised analog digital processor with a microprocessor 
controlled sweep rate and distance windowing function would 
add much to the automation and, consequentially, consistancy 
of results in the TDRL. 

This project is dynamic and holds promise for 
potentially practical applications of applied research. 
There remain many areas to explore by using the OS 
imaginatively. The system is now ready. The true need is 
for the refinement of procedures in using the OS and 
applying it to research areas having merit both instructive 
and practical. 



213 



APPENDIX A 



description op peripheral graphic 3Y3^ ROM 



As Stated in Chapter IV, three peripheral ROM packs have 
been added to the Graphic System to tailor it to the 
specialized needs of the TDRL. Thsy are the Signal 

Processing ROM Pack No. 1, the Signal Processing ROM Pack 
No. 2 (FFT) and the editor. 

The Signal Processing ROM Packs provide fifteen 
mathinatical operations that are essential to advanced signal 
processing and data analysis. Signal processing is 
performed on waveforms or data stored in one-dimensionsal 
arrays. These capabilities do no alter the operation of the 
GS or occupy any of its avaliable RAM space. A summary of 
each ROM command follows. [Ref. 24 25] 

A. Signal Processing ROM Pack No. 1 

1. MAX — finds the maximum of an array. 

2. MIN — finds the minimum of an array. 

3. CROSS -- finis the location of a specified crossing 
level within an array 
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4. DIF2 — performs a two-point dif f erenriation of an 
array . 

5. DIF3 — performs a threa-point dif f arantiation of an 
array. 

6. INT -- integrates an array. 

7. DISP — displays a graph of an array in raw form 
(without axis) . 

3. Signal Processing ROM Pack No. 2 

1. FFT — computes the fast Fouriar transform of an 
array, placing the results in the same array. 

2. IFT -- computes the inverse Fourier transform of an 
array, placing the results in the same array. 

3. CONVL -- convolves two input arrays, placing the 

results in a third array. 

4. CORE -- correlates two input arrays, placing the 

results in a third array. 

5. POLAR — converts an array of FFT data from 

rectangular form (reals and imaginaries) to polar 

form (magnitude and phase) . 

6. TAPER — Multiplies an array by a cosine window of 
program-selectable tapering weights. When tapering 
is selected to be 50%, it provides a Hanning window. 
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7. 



UNLEA V 



sorts an array of interlsaved FFT data into 



two arrays, one containing real and 
imaginary components, 

8. INLSAV — interleaves the real and 
from two input arrays into a third arr 
is acceptable to the IF! command. 

The EDITOR is a ROM device that con 
routines that allow ASCII magnetic files to 
edited off-line, [Ref, 26] 



one containing 

imaginary data 
ay whose format 

tains firmware 
be altered or 
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APPENDIX B 



GENERAL INFORMATION I2S,k O PERA TING SYSTEM 
A. GRAPHING 7ARIABLES 

The following named variables, numeric and constant, are 
used in more than one subprogram. A brief description of 
their purpose and the file in which thay are initialed is 
given below. Note specific program/subroutine for local 
numeric variables, constants and arrays. 



1 1 

VARIABLE 

NAME 


PURPOSE 


1 

FILE 

INITIATED 




STRING CONSTANTS 




E$ 


DATE OF RUN 


2 


?$ 


TARGET DISTANCE 


2 


G$ 


ANTENNA NUMBER 


2 




TARGET NUMBER 


2 


M$ 


VERTICAL SCALE FACTOR - mVOLTS 


2 


?$ 




RUN NUMBER 

— — — ^ — i 


0 

1 



217 



r 

VARIABL2 

NAME 


PURPOSE 


T? T ”7 

INITIATED 


R$ 


"RUN” 


2 


r$ 


HORIZONTAL SCALE FACTOR - SECONDS 


2 


xs 


SPECIAL REMARKS 


2 




NUMERIC CONSTANTS AND VARIABLES 




36 


OPTIMIZATION COMPENSATION 


3 


37 


OPTIMIZATION COMPENSATION 


3 


38 


OPTIMIZATION COMPENSATION 


3 


C5 


TARGET LENGTH 


2 

1 


C6 


TARGET RADIUS 


2 


Cl 


ANTENNA PARAMETER AVERAGING VAL 


2 


L 


OPTIMIZATION FACTOR 


3,4 


L 1 


FLAG 


1 J 

3 


H2 


SIGNALS AVERAGED/WAVEFORM ' 


2 


P7 


FLAG 


2 


?8 

— — — j 


RUN NUMBER 


r 

2 

— 
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r 

VARIABLE 

NAME 


- — , 

PURPOSE 


FILE 

INITIATED 


27 


OPTIMIZATION COMPENSATION 


3 


29 


ENDPOINT DISTANCES OPTICS 


3 


1 

S 1 


MAG, TIME SCALE FACTOR 


2 


S4 


FLAG 


r 

2 


35 


ANTENNA PARAMETER- LENGTH 


3 


36 


ANTENNA PARAM ETER- HI DTH 


3 


1 1 


DPO SHEEP 


2 


V 


MISC CONSTANTS ARRAY 


2,3 

L 


Z1 


r r 

SAMPLING RATE DPO | 2 

_ i 
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APPENDI X C 

LIST OF INPUT PROGMii 

i. GENERAL DESCRIPTION 

1. Location: File 2 

2. Memory space required (less mass storage): 
24552-bytes 

3. Maximum no. of steps: 272 (212 less remarlcs) 

4. Description: INPUT provides the means to initialize 

and drive the TDRL Operating System. It allows input 
of all waveforms, averages acquisitions and stores 
data. It allows for input of antenna target 
parameters and initializes most graphing parameters. 
Inputs are from the keyboard of the 3S, mass storage, 
and File 4. Output is to mass storage and File 3. 
This program is highly operator-processor 
interactive. Numerous messages and directives guide 
the user. 
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B 



LOCAL PBOGRAH ?ARI ABLES AND CONSTANTS 



1. Sianif leant 



— 1 

TYPE 


1 

NAME 


1 

USE 


REMARKS 


ARRAY 


3 


DIWA7 


(512) — TEMPORARY ARRAY FOR 
DIRECT SAVEFORM 


■ ■ 1 
ARRAY 


■ 

C 


AUGWA V 


(512) — TEMPORARY ARRAY FOR 
AUGMENTED SAVSEORM 


I 

ARRAY 


1 

xo 


ACA73 
INWA7 
DIWA7 
ADGNAV 
SCATS AV 


(512) --TEMPORARY ARRAY FOR 
SIGNAL ACQUISITION 
AND AVERAGING 


CONSTANTS 

1 - - - i 


?0 


ALL 


MENU ITEM 



2* Miscellaneous 

a. ARRAY: 30 (512) 

b. STRING CONSTANTS: AS, B$ , C$ , J3, LS , K$ 

Z$, Y$ 

C. NUMERIC VARIABLES AND CONSTANTS: Cl, C2, C3, C4, 

C5 , I, MO, N , N1 , P9, SO, T, X 
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C. MEHOHY SEQOIREMSNTS 



^ • *^05 2 Internal Data Storaae 



r I I 



DESCRIPTIDS 


BYTES 


EACH ACQUIRED WAVEFORM TO BE AVERAGED 


41 14 


AUGMENTED WAVEFORM ARRAY 


4114 


DIRECT WAVEFORM ARRAY 


41 14 


VARIOUS NUMERIC/STRING VARIABLES AND CONSTANTS 


5575 


TOTAL 


17917 


MAXIMUM MEMORY REQUIREMENTS (INTERNAL) FOR INPUT 


42469 1 

i 



2* Ex tern a l M^s Stor age Req uir emen ts (Data Tayje) 



t ~ — 1 

DESCRIPTION 


BYTES 


1 1 

FILE # 


DIRECT WAVEFORM ARRAY 


4114 


50 


INCIDENT WAVEFORM ARRAY 


4114 


51 


AUGMENTED WAVEFORM ARRAY 


4114 

1 


52 

i 


BACKSCATTERED WAVEFORM ARRAY 


4 114 

j 


r ■ ’ " 

1 53 

1 

1 


TOTAL 

_ 1 


r T 

164 561* 
1 



222 



D 



S0B8O0TIHES 



< 1 1 1 T 1 1 



CLASS 


NAME 


LINE 


STEPS 


s 


DESCRIPTION 


DRIVER 


DRIVER 


10 0- 
1380 


129 


40 


INITIALIZE AND DRIVE 
PROGRAM. PASS CONTROL 
TO NEXT FILE 


INPUT 


ACAVG 


2000- 

2730 


74 


23 


ACQUIRES, STORES AND 
AVERAGES DIWAV, INWAV 
AUG WAV INPUTS 


INPUT 


INWAV 


3000- 

3170 


18 
1 9 


6 

6 


INITIALIZES INPUTS. 
STORES INCIDENT WAVE 


INPUT 


DIWAV 


3200- 

3380 


1 9 


6 


INITIALIZES, INPUTS, 
STORES DIRECT WAVE 


INPUT 


AUG WAV 


3400- 

3590 


1 

20 


6 


INITIALIZES, INPUTS, 
STORES AUGMENTED WAVS 


INPUT 
(_ 


SCATWAV 

- 


3600 

3630 


9 


3 


INPUTS TARGET SCATTER 
FROM DIFF. OF INWAV - 
AUGWAV. STORES DATA 
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CLASS 


NAME 


LINE 


STEPS 


% 


1 

DESCRIPTION 1 


SPECIAL 


ANTLOC 


4000- 

4060 


7 


2 


INPOTS ANTENNA PAHA, | 
FROM MASS STORAGE 


SPECIAL 


IGILOC 


4100- 

4160 


7 


2 


TARGET DIMENSIONS 
(IF KNOWN) FROM MASS 
STORAGE 


SPECIAL 


NEHPAR 


4200- 

4330 


4 


1 


INITIALIZES ANTIN AND 
TGTIN FOE STORAGE OF 
NEW ANTENNA AND 
TARGET PARAMETERS 


SPECIAL 


ANIIN 


4400- 

4600 


21 


r 

7 


INPUTS NEW ANTENNA 
PARA. TO MASS STORAGE 


SPECIAL 

L — _J 


TGTIN 


4700- 

4800 




r — 

1 1 


r 

3 

L- - ■ 


INPUTS NEW TARGET 
PARA. TO MASS STORAGE 

.. — _i 
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4430 REN:::::::::::::: driver for NEWPAR to input new antenna 
4440 DIM A5<63> 

4450 A5»0 

4451 FOR I«1 TO €1 STEP 3 
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APPENDIX D 



LIST OF NATH PROGRAM 



A. GENERAL DESCRIPTION 

1. Location: File 3 

2. Memory space required (less mass storage) : 
14544-bytes 

3. Maximum no. of steps: 206 (169 less remarks) 

4. Description: MATH performs the significant 

mathematical computations for the operating system. 
These are fast Fourier and inverse Fourier 
transformations, determination of impulse, srep, and 
ramp responses, noise optimization and compensation, 
and antenna parameter determination. Inputs are from 
File 2, File 4 and mass storage. Output is to File 4 
and mass storage. This program is competely 
automated with no operator-processor interactions. 
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3. LOCAL PROGRAM VARIABLES AND CONSTANTS 
1 . Sign if leant 



1 

TYPE 


1 

NAME 


1 

USE 


' 1 

REMARKS 1 


ARRAY 


30 


OPTIC 

DRIVER 


(512) --METRIC PHYSICAL OPTICS 
FROM RAMP RESPONSE 


ARRAY 


E 


DRIVER 


(257) --FREQUENCY DOMAIN 
BACKSCATTER MAG. 


ARRAY 


EO 


DRIVER 


(257) --FREQUENCY DOMAIN 

IMPULSE RESPONSE PHASE 


ARRAY 


E3 


DRIVER 


(257) --FREQUENCY DOMAIN 
BACKSCATTER PHASE 


ARRAY 


F 


DRIVER 


(257) --FREQUENCY DOMAIN 
DIRECT WAVE MAG. 


ARRAY 


?0 


DRIVER 


(257) —FREQUENCY DOMAIN 

IMPULSE RESPONSE MAG. 


ARRAY 


P3 


DRIVER 


(257) --FREQUENCY DOMAIN 
DIRECT WAVE PHASE 


ARRAY 


lO 


DRIVER 


(512) --TIMS DOMAIN TARGET 
IMPULSE RESPONSE 
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I 

I 



■ 

m 



r 1 

T Y?3 


NAME 


USE 


REMARKS 


ARRAY 


RO 


r ' " '■ ^ "" ' ■ 

DRIVER 


(512) — TIME DOMAIN RAMP 
RESPONSE 


CONSTANTS 


MISC 


ALL 


MISC CONSTANTS AND 
VARIABLES 

. - j 



2 • M iscellane ous 

a. ARRAY: XO (512) ,M (257) ,P (257) 

b. NOMERIC VARIABLES AND CONSTANTS: I, 13, K, M1, M2, 
N, Nl, PI, P2, Q, 31, Q2, Q4, 32, S2, U 

C. HEaORT REQUIREMENTS 

I • iQ.5.2 Internal Memory Require men ts 



I — - 1 — — — I 



DESCRIPTION 


BYTES 


EACH ACQUIRED WAVEFORM TO BE AVERAGED 


4114 


FREQUENCY DOMAIN BACKSCATTER PHASE 


2074 


FREQUENCY DOMAIN BACKSCATTER MAG 


2074 


FREQUENCY DOMAIN DIRECT WAVE PHASE 


2074 


FREQUENCY DOMAIN DIRECT WAVE MAG 


2074 


FREQUENCY DOMAIN IMPULSE PHASE 


2074 
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DESCaiPTIOW 


BYTES 


FREQUENCY DOMAIN IMPULSE MAG 


2074 


EACH ACQUIRED WAVEFORM TO BE AVERAGED 


41 14 


FREQUENCY DOMAIN TEMP STORAGE PHASE 


2074 


FREQUENCY DOMAIN TEMP STORAGE MAG 


2074 


TIME DOMAIN TARGET IMPULSE RESPONSE 


41 1 4 


TIME DOMAIN RAMP RESPONSE 


41 14 


PHYSICAL OPTICS FROM RAMP RESPONSE 


4114 


SCATCH PAD TEMPORARY STORAGE 


4114 


VARIOUS NUMERIC/STRING VARIABLES AND CONSTANTS 


1516 


TOTAL 


34 56 4 


MAXIMUM MEMORY REQUIREMENTS (INTERNAL) FOR INPUT 


49108 

■ j 



D. SXTERNAL MASS STORAGE REQOIREMENTS (DATA TAPE) 






DESCRIPTION 


BYTES 


FILS # 


FREQUENCY DOMAIN IMPULSE RESPONSE 


4 148 


56 


r'" ■ ■■ ' ' . - j i-. . ... 1 

FREQUENCY DOMAIN DIRECT WAVEFORM 


4148 


55 


FREQUENCY DOMAIN BACK SCATTER 


4148 


57 


TIME DOMAIN IMPULSE RESPONSE 


4114 


54 


PHYSICAL OPTICS FROM RAMP RESPONSE 


4 114 


58 


TIME DOMAIN RAMP RESPONSE 


4114 


59 


GRAPHING PARAMETERS 


8364 


60 


TOTAL 


29036 
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SUBROOTI NES 



E. 



— 1 

CLASS 


-I 

NAME 


1 

LINE 


1 

STEPS 


r— 1 

% 


1 — 1 

DESCRIPTION 


DRIVER 


DRIVER 


100- 

1210 


122 


59 


INITIALIZE AND DRIVE 
PROGRAM. PASS CONTROL 
TO NEXT FILE 




FFT 


2000- 

2140 


15 


7 


PERFORMS FAST 
FOURIER TRANSFORM 




OPTIC 


2200- 

2700 


51 


25 


SEPSRATES TARGET 
FROM RANGE NOISE 




ANTDP 


3000- 

3080 


. i. . . . , 1 

9 


4 


UPDATES ANTENNA 
PARAMETERS 




OPTia 


3100- 

3130 


9 

1 


|4 

1 

1 

! 


COMPENSATES PHYSICAL 
OPTICS OPTIMIZATION 

1 
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APPENDIX E 

LIST OF GRAPH PROGRAM 
A. GENERAL DESCRIPTION 



1. 


Location: Fils 


4 




2. 


Memory space 


required (less mass 


storage) : 




44712-bytes 






3. 


Maximum no. of 


steps: 584 (473 less remarks) 


4. 


Descr iption : 


GRAPH provides the means 


to visually 




display the 


results of signal acqui 


sition and 




process ing. 


Subroutines are divided 


into three 




class es; 







a) GENERAL — those subroutines that are applicable to 
most of the other routines in GRAPH. They 
initialize the program, format graphics, title, 
lable, and overlay a grid. 

b) INVERSE S3ATTSRING--thosa routines applicable to 

the inverse scattering solution of time domain 
electromagnetic rransients. They are rhe 

secondary drivers for the development of graphics. 
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c) PRONY’3 TSCHNIQUE-TOUtinss applicable to the 
Prony niethod of solution. 

Inputs to GRAPH are from File 3, File 1 , the GS keyboard 
and mass storage. Outputs are to the GS CRT, mass storage. 
File 2, File 3, File 6, and File 8. This program is highly 
operator-processor interactive. Numerous messages and 
directives guide the user. Signicant graphs are provided 
for operator evaluation of processing results. 

B. LOCAL PSOGRia VARIABLES AND CONSTANTS 
"• • Si gn if ican t 



1 — 1 

TYPE 


NAME 


— 1 

USE 


REMARKS 


ARRAY 


X 


VARIOUS 


(512) — TEMPORARY STORAGE 

AND PROCESSING i 


ARRAY 


Y 


VARIOUS 


(512) — TEMPORARY STORAGE 
AND PROCESSING 


ARRAY 


P 


AXIS 

TICKS 

AXLAB 


(8) —AXIS LABLING 
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' 

TYPE 


NAME 


USE 


1 

REMARKS 1 

1 

j 


CONSTANTS 


SO 


ALL 


FLAG 1 


CONSTANTS 


L$ 


r 

ALL 


HORIZONTAL AXIS LABLE 


CONSTANTS 


NS 


ALL 


VERTICAL AXIS LABLE 


CONSTANTS 


PO 


ALL 


MENU ITEM 


CONSTANTS 


L2 


1 DRIVER 


OPTIMIZATION FACTOR 


CONSTANTS 


AO 


GROP 


HORIZONTAL SCALING 
FACTOR 


CONSTANTS 


A 1 


GROP 


VERTICAL SCALING 
FACTOR 


CONSTANTS 

! 


r* 

W 


GROP 


SCALING VIEWPORT 
COMPENSATOR 

L _ 1 



2 . M isc ellane ous 

a. ARRAY: none 

b. SIRING CONSTANTS: AS, 3$, S$ 

c. NUMERIC VARIABLES AND CONSTANTS: D5 , G1, G2, G3, 

G4, C5, I, 1, Q, PI, P2, ?3, ?4, VI, V2, V3, V4, 
W1, W2, W3, W4, W5, W6 
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C. SEHORT REQOIREMENTS 



1 2 I nt erna 1 Data Stor a 



DESCRIPTION 


1 

BYTES 


VARIOUS TEHPORARI STORAGE/PROCESSING 


3228 


HATRIX LA3LING ARRAY 


82 


VARIOUS NUHSRIC/STRING VARIABLES AND CONSTANTS 


1566 


TOTAL 


9866 


MAXIMUM MEMORY REQUIREMENTS (INTERNAL) FOR INPUT 


54578 



2. Ex ta rn al Hass Sto r aas Requir ement s (Data Tap e) ; no ne 

D. SUBROOTIHES 






CLASS 


NAME 


LINS 


STEPS 1% 

1 

I 1 


DESCRIPTION 


DRIVER 


DRI VER 


100- 

1190 


110 1 1 9 

1 

1 

1 

1 

- . i _ J 


INITIALIZE AND DRIVE 
PROGRAM. PASS CONTROL 
TO NEXT FILE 


GRAPH 


AXIS 


2000- 

2250 


t 1 

26 I 4 

I 

I 

I _ 1_ _ J 


DRIVES AXIS GRAPHING, 


GRAPH 


TICKS 


2300- 

1 


22 Is 
^ 1 , 


AXIS TICKS INTERVALS 

/ 
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CLASS 


NAME 


LINE 


STEPS 


Of 


DESCRIPTION 


GRAPH 


AXLA3 


2700- 

3190 


50 


9 


LA3LE VERTICAL AND 
HORIZONTAL AXIS 


GRAPH 


TITLE 


3200- 

3310 


r "" ' i 

12 


2 


POSITION GRAPH TITLE 


GRAPH 


GRPAR 


3400- 

3870 


48 


8 


PRINTS PARAMETER 
VALUES AS REQUIRED 


GRAPH 


GRID 


3900- 

4010 


r ^ '■ 

12 


2 


DRAWS GRIDS 


INVERSE 

SCATTER 


GRIN 


5000- 

5410 


42 


7 


GRAPHICS FOE DIRECT^ 
INCIDENT, AUGMENTED 
AND 3ACKSCATTER 
TIME DOMAIN WAVES 


INVERSE 

SCATTER 

. 


GRFOUR 


5500- 

6080 


r 

59 


I 

10 


GRAPHICS FOR FREQ 
DOMAIN PHASE 
AND MAGNITUDE OF 
DIRECT, BACKSCATTER 
AND IMPULSE RESPONSE 


INVERSE 

SCATTER 


GROP 


6100- 

5410 


43 


8 

' 


GRAPHICS FOR PHYSICAL 
OPTICS OF TARGET 


INVERSE 

SCATTER 


GRAMP 1 


6600- 1 
6910 

1 1 


1 

32 

i 


5 

j 


GRAPHICS FOR RAMP 
RESPONSE 

1 1 
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CL ASS 


NAME 


LINE 


STEPS 


% 


DESCRIPTION 


PRONY 


PRIMP 


7000- 

7340 


35 


6 


GRAPHICS FOR IMPULSE 
RESPONSE 


PHONY 


PROTHAN 


7400- 

7710 


32 


5 


GRAPHICS FOR PHONY 
TRANSIENT RESPONSE 


PHONY 

1 


PROLES 


7800- 
8300 
1 


51 


9 


GRAPHICS FOR PHONY 
POLES OF IMPULSE 

I 
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APPENDIX F 



IIEICy. 0PSM7ING SYSTEM PROCEDURE 



In this Appendix is presented a typical run in which a 
sphere and a cone are illuminared by the incident 
alecromagnet ic field. The messages presented and the 
results displayed are those as would be encountered during a 
normal run of the system. 

The OS is first used to acquire ani display a normal 
signal (the sphere) . This signal is than optimized to 
remove noise. The results of the effected waveforms are 
then represented. Finally, the system is returned to INPUT 
and a new target, a cone is acquired. The results for the 
cone are presented. Note that the impulse response is 
optimized by a factor of .5 while all other graphs for the 
cone are shown unootimized. 
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